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Abstract 

We compare two different possibilities to include meson-loop corrections in the 
Nambu-Jona-Lasinio model: a strict 1/Ac-expansion in next-to-leading order and a 
non-perturbative scheme corresponding to a one-meson-loop approximation to the 
effective action. Both schemes are consistent with chiral symmetry, in particular 
with the Goldstone theorem and the Gell-Mann-Oakes-Renner relation. The nu- 
merical part at zero temperature focuses on the pion and the p-meson sector. For 
the latter the meson-loop-corrections are crucial in order to include the dominant 
p — > TTTT-decay channel, while the standard Hartree + RPA approximation only 
contains unphysical qq-decay channels. We find that m,r, 7,^5 (^V') a-iid quantities 
related to the p-meson self-energy can be described reasonably with one parameter 
set in the l/A'^c-expansion scheme, whereas we did not succeed to obtain such a fit 
in the non-perturbative scheme. We also investigate the temperature dependence 
of the quark condensate. Here we find consistency with chiral perturbation theory 
to lowest order. Similarities and differences of both schemes are discussed. 
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1 Introduction 



During the last few years one of the principal goals in nuclear physics has been to explore 
the phase structure of QCD. Along with this comes the investigation of hadron properties 
in the vacuum as well as in hot or dense matter. In principle, all properties of strongly 
interacting particles should be derived from QCD. However, at least in the low-energy 
regime, where perturbation theory is not applicable, this is presently limited to a rather 
small number of observables which can be studied on the lattice, while more complex 
processes can either be addressed by chiral perturbation theory or within effective model 
calculations which try to incorporate the relevant degrees of freedom. 

So far the best descriptions of hadronic spectra, decays and scattering processes are 
obtained within phenomenological hadronic models. For instance the pion electromagnetic 
form factor in the time-like region can be reproduced rather well within a simple vector 
dominance model with a dressed p-meson which is constructed by coupling a bare p- 
meson to a two-pion intermediate state |^ . Models of this type have been successfully 
extended to investigate medium modifications of vector mesons and to calculate dilepton 
production rates in hot and dense hadronic matter P]. 

In this situation one might ask how the phenomenologically successful hadronic models 
emerge from the underlying quark structure and the symmetry properties of QCD. Since 
this question cannot be answered at present from first principles it has to be addressed 
within quark models. For light hadrons chiral symmetry and its spontaneous breaking 
in the physical vacuum through instantons plays the decisive role in describing the two- 
point correlators with confinement being much less important. This feature is captured 
by the Nambu-Jona-Lasinio (NJL) model in which the four-fermion interactions can be 
viewed as being induced by instantons. Furthermore the model allows a study of the chiral 
phase transition as well as the examination of the influence of (partial) chiral symmetry 
restoration on the properties of light hadrons. 

The study of hadrons within the NJL model has of course a long history. In fact, 
mesons of various quantum numbers have already been discussed in the original papers 
by Nambu and Jona-Lasinio and by many authors thereafter (for reviews see [|, 0, H)- 

In most of these works quark masses are calculated in mean-field approximation 
(Hartree or Hartree-Fock) while mesons are constructed as correlated quark- ant iquark 
states (RPA). This corresponds to a leading-order approximation in l/Nc-, the inverse 
number of colors. With the appropriate choice of parameters chiral symmetry, which 
is an (approximate) symmetry of the model Lagrangian, is spontaneously broken in the 
vacuum and pions emerge as (nearly) massless Goldstone bosons. While this is clearly 
one of the successes of the model, the description of other mesons is more problematic. 
One reason is the fact that the NJL model does not confine quarks. As a consequence 
a meson can decay into free constituent quarks if its mass is larger than twice the con- 
stituent quark mass m. Hence, for a typical value of m ~ 300 MeV, the p-meson with 
a mass of 770 MeV, for instance, would be unstable against decay into quarks. On the 
other hand the physical decay channel of the p-meson into two pions is not included in 
the standard approximation. Hence, even if a large constituent quark mass is chosen in 
order to suppress the unphysical decays into quarks, one obtains a poor description of the 
p-meson propagator and related observables, like the pion electromagnetic form factor. 
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Similar problems arise if one wants to study the phase structure of strongly interacting 
matter within a mean-field calculation for the NJL model, although this has been done by 
many authors (see e.g. 0, ^, |10|)- In these calculations the thermodynamics is entirely 
driven by unphysical unconfined quarks even at low temperatures and densities, whereas 
the physical degrees of freedom, in particular the pion, are missing. 

This and other reasons have motivated several authors to go beyond the standard 
approximation scheme and to include mesonic fluctuations. In Ref. ||Tl| a quark-antiquark 
p-meson is coupled via a quark triangle to a two-pion state. Also higher-order corrections 
to the quark self-energy [l^ and to the quark condensate 0] have been investigated. 
However, as the most important feature of the NJL model is chiral symmetry, one should 
use an approximation scheme which conserves the symmetry properties, to ensure the 
existence of massless Goldstone bosons. 

A non-perturbative symmetry conserving approximation scheme has been discussed 
in Refs. and []T3|. In Ref. a correction term to the quark self-energy is included in 
the gap equation. The authors find a consistent scheme to describe mesons and show the 
validity of the Goldstone theorem and the Goldberger-Treiman relation in that scheme. 
The authors of Ref. ||15[ use a one-meson-loop approximation to the effective action in a 



bosonized NJL model. The structure of the meson propagators turns out to be the same 
as in the approach of Ref. 
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Based on this scheme various authors have investigated 
the effect of meson-loop corrections on the pion electromagnetic form factor ||16| and on 

and the scalar channel |18[. However, since the numerical 

in these references the exact 



TT-vr scattering in the vector |[I 

evaluation of the multi-loop diagrams is rather involved 
expressions are approximated by low-momentum expansions. 

Another possibility to construct a symmetry conserving approximation scheme beyond 
Hartree approximation and RPA is a strict l/A'^c-expansion up to next-to- leading order. 
Whereas in the approximation scheme mentioned above the gap equation is modified 
in a selfconsistent way, the corrections in the l/A^^c-expansion scheme are perturbative. 
The consistency of the l/A^^c-expansion scheme with chiral symmetry has already been 
shown in Ref. [|l4i . It has been studied in more detail in Refs. |ll9|, E0| . Recently such an 



| . It has been studied in more detail in Refs. 
expansion has been discussed also in the framework of a non-local generalization of the 
NJL model 



21 



In the present paper we compare the results obtained in the non-perturbative scheme 
with those obtained in the 1/iVc-expansion scheme. In vacuum we focus our discussion 
on the pion and the p-meson, calculated with the full momentum dependence of all ex- 
pressions. Within the 1/iVc-expansion scheme the influence of mesonic fluctuations on the 
pion propagator has been examined closely in Ref. 
recent works by Kleinert and Van den Bossche 



19(1 . This was mainly motivated by 



22 1 , who claim that chiral symmetry is 



not spontaneously broken in the NJL model as a result of strong mesonic fluctuations. 
In Ref. we argue that because of the non-renormalizability of the NJL model new 
divergences and hence new cutoff parameters emerge if one includes meson loops. Follow- 
ing Refs. and []T^ we regularize the meson loops by an independent cutoff parameter 
A^f. The results are, of course, strongly dependent on this parameter. Whereas for 
moderate values of Am the pion properties change only quantitatively, strong instabilities 
are encountered for larger values of Am- In Ref. [|19| we suggested that this might be a 
hint for an instability of the spontaneously broken vacuum state. It turns out that the 
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same type of instabilities also emerge in the non-perturbative scheme. This allows for an 
analysis of the vacuum structure and therefore for a more decisive answer to the question 
whether chiral symmetry gets indeed restored due to strong mesonic fluctuations within 
this approximation. 

In any case, in the l/A^^c-expansion scheme the region of parameter values where insta- 



bilities emerge in the pion propagator is far away from a realistic parameter set [20|. We 



used mTT, /tt, (V"^) the p- meson spectral function to fix the parameters. The last one 
is particularly suited, as it cannot be described realistically without taking into account 



pion loops, to fix the parameters. An important result of the analysis in Ref. |]20| was that 
such a fit can be achieved with a constituent quark mass which is large enough such that 
the unphysical gg-threshold opens above the p-meson peak. Since the constituent quark 
mass is not an independent input parameter this was not clear a priori. In this paper we 
will try the same for the selfconsistent scheme. It turns out that it is not possible to find a 
parameter set, where the constituent quark mass comes out large enough to describe the 
properties of the p-meson reasonably. In fact we encounter instabilities in the p-meson 
propagator which are similar to those we found in the pion propagator for large Km- 

The inclusion of meson loop effects should also improve the thermodynamics of the 
model considerably. A first insight on the influence of mesonic fluctuations upon the 
thermodynamics can be obtained via the temperature dependence of the quark conden- 
sate. It has been shown in Ref. that in the selfconsistent scheme the low-temperature 
behavior is dominated by pionic degrees of freedom which is a considerable improvement 
on calculations in Hartree approximation where quarks are the only degrees of freedom. 
Within this scheme the lowest-order chiral perturbation theory result can be reproduced. 
This is also the case for the 1/iVc-expansion scheme which will be demonstrated in the 
last part of this paper. The non-perturbative scheme also allows for an examination of 



the chiral phase transition whereas this is not possible within the 1/iVc-expansion 
scheme. 

The paper is organized as follows. In Sec. ^ we begin with a brief summary of the 
standard approximation scheme used in the NJL model to describe quarks and mesons 
and afterwards present the scheme for describing quantities in next-to-leading order in 
1/N(.. In Sec. ^we discuss the non-perturbative approximation scheme. The consistency of 
these schemes with the Goldstone theorem and with the Gell-Mann Oakes Renner relation 
will be shown in Sec. ^. The numerical results at zero temperature will be presented in 
Sec. ^. The temperature dependence of the quark condensate at non-zero temperature 
within the two above mentioned approximation schemes will be studied in Sec. ^. Finally, 
conclusions are drawn in Sec. ^. 
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2 The NJL model in leading order and next-to-leading 
order in 1/Nc 

2.1 The standard approximation scheme: Hartree + RPA 

We consider the following generalized NJL-model Lagrangian: 

C = - mo)^ + + (V'^75rV')'] - 9. {{i^l^'ri^f + (^^7^75^^^)'] , (2.1) 

where -j/^ is a quark field with N ^ = 2 fiavors and Nc = 3 colors. Qs and (7^ are coupling 
constants with dimension length'^. In contrast to QCD, color is not related to a gauge 
symmetry in this model, but only relates to a counting of degrees of freedom. However, 
if one defines the coupling constants to be of the order l/N^, the large- A^c behavior of the 
model agrees with that of QCD [0, Although we are not interested in the behavior 
of the model for arbitrary numbers of colors in the present article, the 1 /A^^c-expansion is 
introduced for the purpose of book-keeping. This will allow us to take into account mesonic 
fiuctuations in a symmetry conserving way. In order to establish the expansion scheme, 
the number of colors will be formally treated as variable. All numerical calculations will 
be performed, however, with the physical value, A'c = 3. 

In the limit of vanishing current quark masses mo ( "chiral limit" ) the above Lagrangian 
is invariant under global SU{2)l x SU{2)ii transformations. For a sufficiently large scalar 
attraction this symmetry is spontaneously broken. This has mostly been studied within 
the (Bogoliubov-) Hartree approximation.^ 

The Dyson equation for the quark propagator in Hartree approximation is diagram- 
matically shown in Fig. |^. The selfconsistent solution of this equation leads to a momen- 
tum independent quark self-energy Hfj and therefore only gives a correction to the quark 
mass: 

m = mo + T^nim) = mo + 22 '^^9m J Ti [T m S (p) ] . (2.2) 

Usually, m is called the "constituent quark mass". Here S{p) = (^ — m)^^ is the (Hartree) 
quark propagator and "Tr" denotes a trace in color, fiavor and Dirac space. The sum 
runs over all interaction channels M = a, 7r,p, ai with Fo- = 1, F" = Z75r", T'^"' = 
and F^" = 7^75r". The corresponding coupling constants are Qm = Qs for M = cr or 



* Because of the local 4-fermion interaction in the Lagrangian, exchange diagrams can always be cast in 
the form of direct diagrams via a Fierz transformation. This is well known from zero-range interactions in 
nuclear physics. In particular the Hartree-Fock approximation is equivalent to the Hartree approximation 
with appropriately redefined coupling constants. In this sense we call the Hartree approximation the 
"standard approximation" to the NJL model, although in several references a Hartree-Fock approximation 
has been performed. 



-Q 



Figure 1: The Dyson equation for the quark propagator in the Hartree approximation 
(solid lines). The dashed lines denote the hare quark propagator. 
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M = TT and qm = Qv for M = p oi M = ai. Of course, only the scalar channel (M = a) 
contributes in vacuum. One gets 



7^4- ■ (2-3) 

[2TXpp^ — 771^ + le 



In a 1/Nc expansion of the quark self-energy the Hartree approximation corresponds to 
the leading order. Since Qs is of the order l/Nf. the constituent quark mass m, and hence 
the quark propagator are of the order unity. 

For sufficiently large couplings Qs Eq. ( p.3|) allows for a finite constituent quark mass 
m even in the chiral limit. In the mean-field approximation this solution minimizes the 
ground state energy. Because of the related gap in the quark spectrum, one usually refers 
to this equation as the gap equation, in analogy to BCS theory. 

A closely related quantity is the quark condensate, which is generally given by 

m = -^j (^Tr5(p). (2.4) 

In Hartree approximation one immediately gets from the gap equation 

= , (2.5) 

where we have used the superscript (0) to indicate that this corresponds to a Hartree 
approximation. 

Mesons are described via a Bethe-Salpeter equation. Here the leading order in 1/Nc is 
given by a random phase approximation (RPA) without Pauli-exchange diagrams. This 
is diagrammatically shown in Fig. ^. The elementary building blocks of this scheme are 
the quark- ant iquark polarization functions 

nM(g) = -^J 0^,Tr[TM^S{p+l)TM^S{p-l)], (2.6) 

with Tm, M = cr, TT, p, oi as defined above. Again, the trace has to be taken in color, fiavor 
and Dirac space. Iterating the scalar (pseudoscalar) part of the four-fermion interaction 
one obtains for the sigma meson (pion): 

1 - 2gsil^{q) 1 - 2gsil^{q) 

Here a and b are isospin indices and we have used the notation H"''(g) = H^(g) 6ab- 

In the vector channel this can be done in a similar way. Using the transverse structure 
of the polarization loop in the vector channel, 

nM.,-^'(g) = n,(g) T^^ 5ab ; T^^" = {-g^" + ^) , (2.8) 

one obtains for the p-meson 
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Figure 2: The Bethe-Salpeter equation for the meson propagator in the RPA (double line). 
The solid lines indicate quark propagators. 



Analogously, the Oi can be constructed from the transverse part of the axial polarization 
function lia^- As discussed e.g. in Ref. 11^^ also contains a longitudinal part which 
contributes to the pion. Although there is no conceptional problem in including this 
mixing we will neglect it in the present paper in order to keep the structure of the model 
as simple as possible. 

It follows from Eqs. ( |2.6| ) - ( |2.9D that the functions DmIq) are of order l/N^. Their 
explicit forms are given in App. p. For simplicity we will call them "propagators", al- 
though strictly speaking, they should be interpreted as the product of a renormalized 
meson propagator with a squared quark- meson coupling constant. The latter is given by 
the inverse residue of the function Z)m(q'), while the pole position determines the meson 
mass: 

n-l/ M n -2(0) dllMiq) . /oin\ 

^A/(9)L2 2(0) = 0, g'' = ,0 L2 2(0) . (2.10) 



Again the superscript (0) indicates that m^^^-* and g^f^^^ are quantities in RPA. One easily 
verifies that they are of order unity and 1 / respectively. 



2.2 Next-to-leading order corrections 

With the help of the gap equation, Eq. (pl3|), one can show that the "standard scheme", 
i.e. Hartree approximation + RPA, is consistent with chiral symmetry. For instance, 
in the chiral limit pions are massless, as required by the Goldstone theorem. Of course 
one would like to preserve this feature when one goes beyond the standard scheme. One 
way to accomplish this is to perform a strict 1/A^c expansion, systematically including 
higher-order corrections. In this subsection we want to construct the quark self-energy 
and the mesonic polarization functions in next-to-leading order in 1/A^c- 
The correction terms to the quark self-energy, 

SE{p) = ^S^'^) + 5S(^)(p) , (2.11) 

are shown in Fig. |^. In these diagrams the single lines and the double lines correspond to 
quark propagators in the Hartree approximation (order unity) and to meson propagators 
in the RPA (order 1/Nc), respectively. Recalling that one obtains a factor A''^ for a closed 
quark loop one finds that both diagrams are of order the l/N^. One can also easily 
convince oneself that there are no other self-energy diagrams of that order. 

According to Eq. (|2.4| ), the l/A^^c-correction to the quark condensate is given by 

6{^^) = -zj (|^Tr55(p), (2.12) 



7 




Figure 3: The 1 / Nc- corrections ^S*^"^ (^^ft) o,nd (right) to the quark self-energy. 

with 

5S{p) = Sip)5i:ip)Sip) (2.13) 

being the 1/A^c-correction to the Hartree quark propagator S{p). Since we are interested 
in a strict 1/Nc expansion, the self-energy correction must not be iterated. 
The l/A'c-corrected mesonic polarization diagrams read 

IlMiq) = HaM) + ^ni?(g) . (2.14) 

k=a,b,c,d 

The four correction terms SH^^j to SH^^ together with the leading-order term Um are 
shown in Fig. ^. Again, the lines in this figure correspond to Hartree quarks and RPA 
mesons. Since the correction terms consist of either one RPA propagator and one quark 
loop or two RPA propagators and two quark loops they are of the order unity, whereas 
the leading-order term is of the order Nc- 

In analogy to Eqs. ( |2.7D , (|2.9| ) and ( |2.10D the corrected meson propagators are given 



by 

DM{q) = ^ ^ . . , (2.15) 

while the corrected meson masses are defined by the pole positions of the propagators: 

D^Kq)\,2=ml, = 0. (2.16) 



As we will see in Sec. [4.2| this scheme is consistent with the Goldstone theorem, i.e. in the 



chiral limit it leads to massless pions. Note, however, that because of its implicit definition 
mjvf contains terms of arbitrary orders in l/N^, although we start from a strict expansion 
of the inverse meson propagator up to next-to-leading order. This will be important in 
the context of the Gell-Mann-Oakes-Renner relation. 

For a more explicit evaluation of the correction terms it is advantageous to introduce 
the quark triangle and box diagrams which are shown in Fig. |^. The triangle diagrams 
entering into ^n^-* and SH^^J can be interpreted as effective three-meson vertices. 

For external mesons Mi, M2 and M3 they are given by 

^^-y^jTr [TMjS{k)TM2iS{k -p)TM3iS{k + q)] 

+Tr [^M,^Sik - q)^M,^S{k+p)^M,^S{k)]] (2.17) 



with the operators Tm as defined below Eq. (|2.2| ). We have summed over both possible 
orientations of the quark loop. For later convenience we also define the constant 
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O <o> 




RPA 



(a) 



(b) 



(d) 



Figure 4: Contributions to the mesonic polarization function in leading (RPA) and next- 
to-leading order in l/Nc- 



1 f d'^p 



(2.18) 



which corresponds to a quark triangle coupled to an external scalar vertex and a closed 
meson loop. 

The quark box diagrams are effective four-meson vertices and are needed for the eval- 
uation of ^n^^ and SIV^j . If one again sums over both orientations of the quark loop they 
are given by 



— ^Fmi ,M2 ,M3 ,M4 (Pl , P2 , ) 
d^k 



Tr [TM^iS{k)TM2iS{k - p2)TMjS{k - p2 - P3)TMjS{k + pi)] 
-Tr [TMjS{k ~ pi)TMjS{k + p2 + p3)TM3'>'S{k p2)TM2iS{k)]) (2.19) 



With these definitions the various diagrams can be written in a relatively compact 
form. For the momentum independent correction term to the quark self-energy we get 



^^.(0) I ^,Dm{p)Tm,mAp,-p) = I^.(0)A 



(2.20) 



In principle there should be also a sum over the quantum numbers of the meson which 
connects the quark loop with the external quark legs, but all contributions from other 
mesons than the cr- meson vanish. The factor of 1/2 is a symmetry factor which is needed 
because otherwise the sum over the two orientations of the quark propagators, which is 
contained in the definition of the quark triangle vertex (Eq. (|2.17|) ) would lead to double 
counting. 



P3,M3 

p' = -p-q,Mz /< Pi , ^\ Pi , Mi 

Figure 5: (Left) The quark triangle vertex —iV Mi,M2,M3iliP)- (Right) The quark box vertex 

— '^^Ml,M2,M3,M4(Pl)P2,P3)■ 
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The evaluation of the momentum dependent correction term ST,^''^ is straight forward: 
5EW(A;) = [ -^,DM{p)TMS{k~p)TM. (2.21) 

Inserting these expressions for 5Tj^°^ and ^S*^*-' into Eq. ( |2.13| ), the l/A^^c-correction term 
to the quark condensate, Eq. ( p.l2| ), can be brought into the form 

,«.^.) . -2M^ . (2.22) 
For the mesonic polarization diagrams we get: 

^n^l?) = ^ /" 7^4i ^ rM,Mi,M2(?,p)^A/i(p)rM,Mi,M2(-g, -p)^M2(-p-g) , 

J (2vr)4 ^ 

^n£^(?) = ^ rM,M,a(g, -g) ^a(o) /" 7^ ^ '^mi,miAp^-p) DmAp) , 

J y ^1 Ml 

= -zrA/,Af,a(g,-g)/^a(0)A . (2.23) 

The symmetry factor of 1/2 for 511 and ^n^"^^ has the same origin as in Eq. (|2.20|) . 
Similarly in 511^"^ we had to correct for the fact that the exchange of Mi and M2 leads to 
identical diagrams. 

For the further evaluation of Eqs. (|2.2CI|) to (|2.23|) we proceed in two steps. In the 



first step we calculate the intermediate RPA meson-propagators. Simultaneously we can 
calculate the quark triangles and box diagrams. One is then left with a meson loop which 
has to be evaluated in a second step. 

The various sums in Eqs. ( p.20| ) to ( |2.23| ) are, in principle, over all quantum numbers 



of the intermediate mesons. However, for most applications we expect that the most 
important contributions come from the pion, which is the lightest particle in the game. 
For instance, the change of the quark condensate at low temperatures should be dominated 
by thermally excited pions. Also, for a proper description of the p-meson width in vacuum 
we only need the two-pion intermediate state in diagram SU^^j . Other contributions to this 
diagram, i.e vrai, pa, pp and aiOi intermediate states, are much less important since the 
corresponding decay channels open far above the p-meson mass and - in the NJL model 
- also above the unphysical two-quark threshold. Hence, from a purely phenomenological 



point of view, it should be sufficient for many applications to restrict the sums in Eq. ( |2.23|) 
to intermediate pions. However, in order to stay consistent with chiral symmetry, we 
have to include intermediate sigma mesons as well. On the other hand, vector- and axial- 
vector mesons can be neglected without violating chiral symmetry. Since this leads to 
an appreciable simplification of the numerics we have restricted the intermediate degrees 
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of freedom to scalar and pseudoscalar mesons in the present paper. Of course, in order 
to describe a p-meson, we have to take vector couphngs at the external vertices of the 
diagrams shown in Fig. ^. 

3 Non-perturbative symmetry conserving schemes 
3.1 Axial Ward identities 

The disadvantage of the l/Ai'c-expansion scheme is that it is perturbative. Although 
we have constructed the 1/iVc corrections to the Hartree quark self-energy (Fig. ^ we 
did not selfconsistently include such diagrams in the gap equation. Since the iteration 
would produce terms of arbitrary orders in one is not allowed to do so in a strict 

expansion scheme. Therefore, all correction diagrams we have discussed in the previous 
section consist of "Hartree" quark propagators. This perturbative treatment should work 
rather well as long as the l/iV^ corrections to the quark self-energy are small compared 
with the leading order, i.e. the constituent quark mass. On the other hand it is clear that 
the scheme must fail to describe the chiral phase transition, e.g. at finite temperatures. 
Here a non-perturbative treatment is mandatory. 

Therefore, in this section, we want to follow a different strategy, exploiting the fact 
that the Goldstone theorem is basically a consequence of Ward identities: Consider an 
external axial current j|^5 coupled to a quark. Then, in the chiral limit, the corresponding 
vertex function F^g is related to the quark propagator S{p) via the axial Ward-Takahashi 
identity 

?'^r«5(p,g) = S'\p + q)^,T^ + l,r^S~\p), (3.1) 

where p and p + q are the 4-momenta of the incoming and outgoing quark, respectively. 
Obviously, for a non- vanishing constituent quark mass, the r.h.s. of this equation remains 
finite even for g — 0. Consequently F°5(p, g) must have a pole in this limit, which can 
be identified with the Goldstone boson. Moreover the explicit structure of the Goldstone 
boson can be constructed from the structure of the axial vertex function. 

As a first example, let us start again from the Hartree gap equation (Eq. (|2.2| ), Fig. |l|) 
and construct the axial vertex function by coupling the propagator to an external axial 
current. This is illustrated in Fig. ^ In the upper line, the first term on the r.h.s. 
describes the coupling to the bare quark, corresponding to a bare vertex 7^75r". In the 
second term, however, the current is coupled to a dressed quark, and therefore we have 
to use the same vertex function as on the l.h.s. of the equation: 

r^sb,?) = l.l^r^ + Y.'^^QmTm j -0y,Tv[TMS{k + q)T%{k,q)S{k)]. (3.2) 

Here S{k) denotes the quark propagator in the Hartree approximation. As in Eq. ( |2.2| ) 
the sum runs over all interaction channels, but of course only the pseudoscalar and the 
axial vector contributions do not vanish. Contracting Eq. ( p.2[ ) with q^ one obtains a 
linear equation for g^F^g. One can easily verify that in the chiral limit the solution of 
this equation is given by the axial Ward-Takahashi identity, Eq. (|3.1|) . To this end we 
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Figure 6: Vertex function for an external axial current coupled to a "Hartree" quark. 

replace F^g on both sides of the equation by the expressions given by Eq. (|3.1| ) and 
check whether the results agree. On the r.h.s. the insertion of Eq. basically amounts 
to removing one of the quark propagators from the loop. In this way the loop receives the 
structure of the quark self-energy and we can use the gap equation, Eq. ( |2.2| ) to simplify 
the expression. For mo = the result turns out to be equal to the l.h.s. of the equation, 
which proves the validity of the axial Ward-Takahashi in this scheme. 

We have seen above, that this implies the existence of a massless Goldstone boson in 
the chiral limit. As illustrated in the second and third line of Fig. |], the selfconsistent 
structure of Eq. (|3.2| ) for the dressed vertex F^g leads to an iteration of the quark loop 
and an RPA pion emerges. Hence we can identify the Goldstone boson with an RPA pion. 

Obviously the above procedure can be generalized to other cases: Starting from any 
given gap equation for the quark propagator we construct the vertex function to an 
external axial current by coupling the current in all possible ways to the r.h.s. of the 
equation. As long as the gap equation does not violate chiral symmetry this automatically 
guarantees the validity of the axial Ward-Takahashi identity and therefore the existence 
of a massless pion in the chiral limit. The structure of this pion can then be obtained 
from the structure of the vertex correction. 

As an example we start from the extended gap equation depicted in the upper part 
of Fig. 1^. There, in addition to the Hartree term, the quark is dressed by RPA mesons. 
These are defined in the same way as before (Fig. |^), but now selfconsistently using the 
quark propagator which results from the extended gap equation. Therefore the RPA pions 
are no longer massless in the chiral limit. However, following the strategy described above 
we can construct the consistent pion propagator. To that end we couple again an external 
axial current to both sides of the gap equation. The resulting equation for the vertex 
function is also shown in Fig. |^ (middle part). The additional term in the gap equation 
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leads to two new diagrams which were not present in Fig. ^: In the first the current 
couples to a quark-antiquark loop of the RPA meson while in the second it couples to 
the quark inside of the meson loop. Again, one can easily check that the vertex function 
and the quark propagator fulfill the axial Ward-Takahashi identity Eq. ( p.l|) in the chiral 
limit. 

In principle one can construct the corresponding massless Goldstone boson from the 
quark-antiquark T-matrix given in the lower part of Fig. 0. In practice, however, this 

- . -Q_ . -(Du 




Figure 7: Self consistent scheme with a non-local self-energy term: Gap equation (upper 
part), equation for the vertex function of an external current (middle) and the correspond- 
ing equation for the quark-antiquark T-Matrix (lower part). The double line denotes an 
RPA meson propagator (see Fig. which is self consistently constructed from the dressed 
quark propagators of the present equation (solid line). 
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is very difficult. In fact, already the solution of the extended gap equation is difficult, 
since the additional self-energy term is non-local, leading to a non-trivial 4-momentum 
dependence of the quark propagator. Note that this propagator has to be selfconsistently 
used for the calculation of the RPA-meson propagator. Therefore the authors of Ref. 
suggested to drop the non-local terms, but to keep a particular class of local diagrams 
which arises from the combined iteration of the quark loop and the meson loop. This gap 
equation is shown in Fig. |. Because of the restriction to local self-energy insertions, we 
will call this scheme the "local selfconsistent scheme" (LSS). It will be discussed in the 
next subsection. 



3.2 The local selfconsistent scheme 

The gap equation for the constituent quark mass in the local selfconsistent scheme (upper 
part of Fig. H) reads 

m = rjiQ + S(m) = mo + S//(m) + 5S(m) . (3.3) 

Here T,h is the Hartree contribution to the self-energy as defined in Eq. ( p.2|) . The 
correction term 5S corresponds to the third diagram on the r.h.s. of Fig. p. We have 
explicitly indicated that the self-energy diagrams have to be evaluated selfconsistently at 
the quark mass m, which comes out of the equation. Because of the new diagram 5S, this 
mass is in general different from the Hartree mass. However, since all diagrams in the 
LSS are constructed from the constituent quarks of Eq. (p^), we prefer not to introduce a 

- -Q .a_ 




Figure 8: The "local selfconsistent scheme": Gap equation (upper part), equation for the 
vertex function of an external current (middle) and the corresponding equation for the 
consistent meson propagator (lower part). The double line denotes an RPA meson propa- 
gator (see Fig. which is selfconsistently constructed from the dressed quark propagators 
of the present equation (solid line). 
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new symbol for this mass. This has the advantage that we can also keep the notations for 
the quark propagator S{p) = {jf — m)^^ , quark-antiquark loops, triangles etc. which we 
introduced earlier. The general structure of these diagrams is the same in all schemes we 
discuss in this article. Therefore we introduce the convention that in the l/A'^c-expansion 
scheme m denotes the Hartree mass, while it denotes the solution of Eq. ( p.3| ) in the LSS 
and all diagrams should be evaluated at that mass, unless stated otherwise. 

The self-energy term 5S consists of a quark loop dressed by an RPA-meson loop. The 
quark loop is coupled to the external quark propagators via the NJL point interaction. It 
can again be shown, that only the scalar interaction contributes. Hence (5S is given by 

5t = -2gsA, (3.4) 

where A is the constant defined in Eq. ( |2.18| ). 

Because of this additional self-energy diagram in the gap equation, the RPA is not 
the consistent scheme to describe mesons: In the chiral limit, RPA pions are no longer 
massless. Hence, in order to find the consistent meson propagators, we proceed in the 
way discussed in the previous subsection. 

The equation for the axial vertex function is shown in the middle part of Fig. |^. Com- 
pared with the corresponding equation which follows from the Hartree approximation 
(Fig. ^) there are three extra terms. This leads to three additional polarization dia- 
grams, which have to be iterated in the Bethe-Salpeter equation for the consistent meson 
propagator (lower part of Fig. |). 

Obviously these diagrams are identical to SU^^j , ^-^id SU^^j , which we defined in 
Sec. (Fig. Eq. ( |2.23|) ), i.e. the new meson propagators are given by 

DM{q) = ^ . . , (3.5) 

with 

riM(g) = nM(g) + ^n^/ (^) • (3-6) 

This structure agrees with the result of Ref. W^- In that reference the scheme was 



motivated by a l/Ai'c expansion. However, one should stress again, that the selfconsistent 
solution of the gap equation mixes all orders in XjNc- Moreover, the next-to-leading order 
self-energy correction term ST^^^ (cf. Fig. ^ is not contained in the gap equation of Fig. H. 
Therefore the consistency of the scheme cannot be explained by l/A'^c-arguments. In fact, 
our discussion shows that the structure of the consistent pion propagator can be derived 
from the gap equation without any reference to l/Ai'c counting. 

On the other hand, if one performs a strict l/iV^ expansion of the mesonic polarization 
diagrams up to next-to-leading order one exactly recovers the diagrams shown in Fig. ^ 
||14|| . This is quite obvious for the diagrams 5H^ to 5Hjy'^\ which are explicitly contained 



in Eq. ( p.6| ). Diagram 5hS5\ which seems to be missing, is implicitly contained in the 



'-M ' 

quark-antiquark loop via the next-to-leading order terms in the quark propagator, which 
arise from the extended gap equation. 

In this sense, the LSS may be viewed as the simplest non-perturbative extension of 
the standard scheme which is consistent with the Goldstone theorem and which contains 
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all mesonic polarization diagrams up to next-to- leading order in l/N^. However, since 
the diagram ^E^'') is not contained in the gap equation, this is not true for the quark 



condensate: If we evaluate Eq. (2^) with the quark propagator of the present scheme, we 
obtain 

m = = -^^^ - A . (3.7) 

Performing a strict l/Nf. expansion of this expression and only keeping the next-to-leading 
order term one does not recover Eq. ( p.22|) but only the contribution of 6'E^"-\ This might 



be the reason why the authors of Ref. WM determine the quark condensate as 



/7/\ ^ m-mo 

2gs 2gs 

In contrast to Eq. (|3.7|) this expression reduces to the perturbative result in a strict 
expansion. Moreover, as we will discuss in Sec. [4.2| , it is consistent with the Gell-Mann- 
Oakes-Renner relation. On the other hand, Eq. (|3.8| ) does obviously not follow from 



c 



Eq. (|2.4|) with the quark propagator of the present scheme. A possible resolution to this 



problem was given in Ref. [|I5[, where the LSS was derived using functional methods. The 
meson propagators obtained in that way are identical to Eqs. ( |3.5| ) and (|3.6|), while the 
quark condensate is given by Eq. ( |3.8|) . This will be briefly discussed in the following 
subsection. 

Finally, we would like to comment on the name "local selfconsistent scheme" , which we 
have introduced in order to distinguish this scheme from the perturbative l/N^ expansion. 
We call this scheme "selfconsistent" because the quark propagator, which is determined 
by the gap equation is selfconsistently used in the loops and the RPA-meson propagator 
on the r.h.s. of that equation. However, as we have seen, the scheme is not selfconsistent 
with respect to the mesons: The improved meson propagators given by Eqs. (|3.5| ) and 
( ^.6[ ) are different from the RPA mesons which are used in the gap equation and hence 
as intermediate states in the mesonic polarization functions 5S^") to 5Tj^^\ On the other 
hand, if we had used the improved mesons already in the gap equation, our method of 



Sec. would have led to further mesonic polarization diagrams in order to be consistent 
with chiral symmetry. Obviously, the construction of an expansion scheme, which is 
selfconsistent for quarks and mesons, is an extremely difficult task. 

3.3 One-meson-loop order in the effective action formalism 

Both, the non-local selfconsistent scheme, which we briefly discussed in Sec. ^]T] (Fig. |^, 
and the local selfconsistent scheme can be derived from functional methods: The non- 
local selfconsistent scheme can be obtained as a ^-derivable theory p5| , p6| if one includes 



the "ring sum" in the generating functional. The present section is devoted to a brief 
discussion on how the local selfconsistent scheme can be derived from a one-meson-loop 
approximation to the effective action. The interested reader is referred to Refs. |27[| . 
Here we will basically follow Ref. ||15|| . 

In this section we drop the vector and axial vector interaction and start from a La- 
grangian which contains only scalar and pseudoscalar interaction terms: 

C = ^{t^ - mo)^ + gs [{^^'f + {^ti5T^y] . (3.9) 
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The partition function of the system can be expressed in terms of the path integral 

Z = e-^ = J ^ (3,10) 

with the Euchdean action 
Ii^P\^|J) = [ d^XE |^S(5r7o-^7■ V + mo)V'-^?s((^S^)' + (V'S^75r^)')| • (3.11) 



The integration is here over a Euchdean space-time volume d^XE, where dr corresponds 
to idt- The standard procedure is now to bosonize the action by introducing auxiliary 
fields $l,a = {0,1,2,3} 



V{^^)V{^)Vm exp I - ^) - i7 y d'xEiK + 2gs^hora^y\ , (3.12) 

with Ta = (Iji'-f^r). Then the action contains only bilinear terms in the quark fields, 
so that they can be integrated out. After performing a shift of the auxiliary fields, 
$a = $^ + (mo,0), one finally arrives at the bosonized action 

/($) = -Tr In + / rfV(*^ -2mo$o + m2) , (3.13) 
where is the Dirac operator 

S-'=^0dr-ll-V + Ta^a. (3.14) 

The symbol Tr in Eq. ( p.l3| ) is to be understood as a functional trace and a trace over 
internal degrees of freedom like fiavor, color and spin. Tr In is the quark-loop contri- 
bution. The imaginary part of this term vanishes for the SU (2) case and we can rewrite 
the action as: 

/($) = ~^TrlnS-^^S-^ + 7^ / d^^Ei^^ - 2mo$o + ml) . (3.15) 

The effective action r($) is defined as Legendre transform of the generating functional 
W{i). Its stationary point ($a), i.e. 



5r($) 



5$n 



= 0, (3.16) 



represents the vacuum expectation values of the fields. 

The quark condensate can be expressed via the expectation value of $o. It is given as 

dW 1 

«'« = a^ = -25:«*°>-'"°)- 

Another important feature of the effective action is that the inverse propagators of the 
fields (in our case the propagators for vr- and cr-mesons) can be generated in a symmetry 
conserving way by second-order derivatives 
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To obtain an expression for the effective action the path integral is evaluated using 
the saddle point approximation. The lowest-order contribution to the effective action is 

r(o)($) = /($) . (3.19) 

This corresponds to the mean-field (Hartree) approximation |2^. The vacuum expectation 
values of the fields in mean-field approximation coincide with the stationary point of the 
action /($). This is obvious if one combines Eq. ( p. 161 ) and Eq. ( |3.19| ). Including quadratic 



mesonic fluctuations leads to the following expression for the effective action |27| : 



r(*) = /(*) + . (3.20) 

The second term in the above expression contains the mesonic fluctuations. As discussed 



in Ref. the method is only meaningful if the second-order functional derivative which 
enters into this term is positive definite. Otherwise severe problems arise due to an ill- 
defined logarithm, which would then be complex. We will come back to this point in 



Sec. 5.2 



Determining the stationary point of the effective action in Eq. (|3.20D leads to the 



following "gap equation" 



($0) - mo - ^hH^^o)) - St{{%)) = . (3.21) 



Here T,h and 5S are the same functions we already defined in Eqs. ( |2.2| ) and ( |3.4| ) in the 
context of the Hartree- and the LSS gap equation. In fact, Eq. ( |3.21|) is identical to the 
LSS gap equation, Eq. (|3.3|), if we identify ($0) with the LSS-constituent quark mass m. 

In the same way we exactly recover the meson structure of the LSS if we evaluate 
Eq. ( 3.181 ) at the stationary point. This means, the "local selfconsistent scheme" which 



was constructed from a somewhat arbitrary starting point in Sec. p.2| can be derived in a 
systematic way in the effective action formalism. However, the interpretation is different: 
As emphasized in Ref. |]15| , the solution of the gap equation is only the expectation value 



of the $0 field and does not correspond to the pole of the quark propagator. This becomes 
clear if we look at the quark condensate, which is given by Eq. ( |3.17D . The r.h.s. of this 



equation is identical to Eq. (|3.8|) and therefore different from Eq. ( p. 71) , which was derived 
by taking the trace over what we called the "quark propagator" in Sec. P^. 



Hence, within the effective action formalism, Eq. ( |3.8|) is the correct expression for 
the quark condensate (in that approximation scheme), whereas the gap equation should 
not be interpreted as an equation for the corresponding inverse quark propagator. In the 
following, we will adopt this point of view. For simplicity, however, we will still call m 
a "constituent quark mass" and (^ — m)~^ a "quark propagator", although this is not 
entirely correct. 

4 Consistency with chiral symmetry 

By construction, the LSS is consistent with axial Ward-Takahashi identities and hence - as 



discussed in Sec. - with the Goldstone theorem. Since the mesonic polarization func- 



tions of the LSS contain all diagrams up to next-to-leading order of the l/A^^c-expansion 



18 



scheme and the various contributions to the pion mass have to cancel order by order in 
the chiral hmit, this imphes that the l/N^ scheme discussed in Sec. p.2| is also consistent 
with the Goldstone theorem. 

Nevertheless, for the numerical implementation it is instructive, to show the consis- 
tency of the different schemes with chiral symmetry on a less formal level. Since most of 
the integrals which have to be evaluated are divergent and must be regularized one has 
to ensure that the various symmetry relations are not destroyed by the regularization. 
To this end, it is important to know how these relations emerge in detail. This will also 
enable us to perform approximations without violating chiral symmetry. As we will see 
in Sec. |5.2| , this is very important for practical calculations within the LSS, which cannot 
be applied as it stands. 

For both, the l/A'c-expansion and the LSS, we begin our discussion with the explicit 
proof of the Goldstone theorem. This was given first by Dmitrasinovic et al. fl^. Af- 



ter that we discuss the Gell-Mann-Oakes-Renner (GOR) relation. This is of particular 
interest in the context of the proper definition of the quark condensate in the LSS (cf. 
Eqs. Q and i^). 



4.1 1/A^c-expansion 

We begin with the 1/iVc-expansion scheme. For the Goldstone theorem one has to show 
that, in the chiral limit, the inverse pion propagator vanishes at zero momentum, 

2c/,n^(0) = 1 for mo = 0. (4.1) 

As before we use the notation 11^^ = 5ah^-K- The function 11^* has been defined in 
Eq. (|2.14|) . It consists of the RPA polarization loop 11"^ and the four l/A^^c-correction 
diagrams ^Ili*'''"*, k = a,b, c, d. Restricting the calculation to the chiral limit and to zero 
momentum simplifies the expressions considerably and Eq. ( |4.1| ) can be proven analyti- 
cally. 

For the RPA loop one obtains 

2g,U^{0) = ^. (4.2) 
m 

This is the relation which guarantees the consistency of the Hartree -|- RPA scheme: 
In Hartree approximation we have m = tuq + Hh and hence Eq. ( |4.1|) , is fulfilled by 
Eq. (|4.2| ). Since the gap equation is not changed in the perturbative l/N^ expansion, this 
remains true, if we include the next-to-leading order. Therefore we have to show that the 
contributions of the correction terms add up to zero: 

= for mo = 0. (4.3) 

k=a,b,c,d 

The correction terms 6111 are defined in Eq. ( p.23| ). Let us begin with diagram 6Ui^\ As 
mentioned above, we neglect the p and ai subspace for intermediate mesons. Then one can 
easily see that the external pion can only couple to a na intermediate state. Evaluating 
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the trace in Eq. ( p.l7| ) for zero external momentum one gets for the corresponding triangle 
diagram: 

K'^,A0,P) = -Sab ^N,Nf 2m I{p) , (4.4) 
with a and b being isospin indices and the elementary integral 

^^^^ " / (27r)4 - m2 + te){{k + - m? + le) ' ^^'^^ 
Inserting this into Eq. ( p.23| ) we find 

SU^:^'^\0) = i5ab I ^^{m,NfI{p))Hm' D^ip) D^ip) . (4.6) 

Now the essential step is to realize that the product of the RPA sigma- and pion propa- 
gators can be converted into a difference [n], 

to finally obtain 

^nwa^(o) = -s^, m^Nf I ^,^np){DAp) - dAp)} ■ (4.8) 

The next two diagrams can be evaluated straightforwardly. One finds: 

5nWa6(0) = -Sat I ^,{D<rip) (Ap) + ^(0) " (/ - ^m') K{p)) 

+D^{p) (Slip) + 3/(0) - 3/ Kip)) } , 
Sui^)'^\0) = -SabAN^Nf j ^J{p){-D^{p)-D^{p)]. (4.9) 

The elementary integral K{p) is of the same type as the integral I{p) and is defined in 
App. 0. 

Finally we have to calculate ^Ili- (0). According to Eq. ( p.23| ), it can be written in 
the form 

^nWa^'(o) = -iV%^^{Q,Q)D,{Q)/\. (4.10) 
For the constant A, defined in Eq. (|2.18|), one obtains 



/d 
^{ D^{p) (2 I{p) + /(O) - (/ - 4m^) K{p)) 

+D^{p){?,m - Sp'Kip)) } . (4.11) 



Evaluating -D(j(0) in the chiral limit and comparing the result with Eq. ( [4. 4] ) one finds 
that the product of the first two factors in Eq. ( [4.10| ) is simply 6ab/m, i.e. one gets 

5Ui^)-\0)=S^^^ . (4.12) 
m 

20 



With these results it can be easily checked that Eq. ( [4 .31 ) indeed holds in this scheme. 

As already pointed out, most of the integrals we have to deal with are divergent and 
have to be regularized. Therefore one has to make sure that all steps which lead to 
Eq. (|4.3| ) remain valid in the regularized model. One important observation is that the 
cancellations occur already on the level of the p-integrand, i.e. before performing the 
meson-loop integral. This means that there is no restriction on the regularization of this 
loop. We also do not need to perform the various quark loop integrals explicitly but we 
have to make use of several relations between them. For instance, in order to arrive at 
Eq. ( [4.12|) we need the similar structure of the quark triangle r^_7r^o-(0, 0) and the inverse 



RPA propagator Df^{0)~^. Therefore all quark loops, i.e. RPA polarizations, triangles 
and box diagrams should be consistently regularized within the same scheme, whereas 
the meson loops can be regularized independently. 

Going away from the chiral limit the pion recieves a finite mass. To lowest order in 
the current quark mass it is given by the Gell-Mann-Oakes-Renner (GOR) relation, 

mlf^ = -mo{^^) . (4.13) 

However, in the l/A'^c-expansion scheme we cannot expect, that the GOR relation holds 
in this form. In Sec. ^ we have calculated the quark condensate in leading order and 
next-to-leading order in l/Nc- Hence, to be consistent, we should also expand the l.h.s. 
of the GOR relation up to next-to-leading order in l/Nc- 

mT^ff^ + ^T^fl + 5^lfT^ = -^0 ( m^'^ + m)) ■ (4.14) 



Here, similar to the notations we already introduced for the quark condensate, rrt^^^ and 
f'^^^^ denote the leading order and 5m^ and 5/^ the next-to-leading order contributions to 
the squared pion mass and the squared pion decay constant, respectively. Since the GOR 
relation holds only in lowest order in mo, Eq. (^4. 14|) corresponds to a double expansion: 
m^ has to be calculated in linear order in mo, and {^ip) in the chiral limit. 

The leading-order and next-to-leading-order expressions for the quark condensate are 
given in Eqs. ( p^.5| ) and ( 2.22 ). The pion decay constant is calculated from the one- 



pion to vacuum axial vector matrix element. Basically this corresponds to evaluating the 
mesonic polarization diagrams. Fig. ^ coupled to an external axial current and to a pion. 
This leads to expressions similar to Eqs. (|2.6| ) and ( |2.23|) , but with one external vertex 



equal to 7'^75^, corresponding to the axial current, and the second external vertex equal 
to g-Kqqi'^^T^ 1 corresponding to the pion. Here the 1/A^c-corrected pion-quark coupling 
constant is defined as 



9-1 = + ^9-1 = ^^U^-mi , (4.15) 



analogously to Eq. ( |2.10| ). Now we take the divergence of the axial current and then use 
the relation 

751^ = 2m75 +-i^S-\k + p) +S-\k)-f^ (4.16) 
to simplify the expressions ||l4l. One finds: 



( n^(g) - n^(o) , - n.(o) ^ ^ ^ 

A = 9nqq m + /^a(0) A ) ^ ^ . (4.17) 
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In the chiral limit, 



m: 



0, Eqs. ( |2.10| ) and ( |4.15| ) can be employed to replace the 



difference ratios on the r.h.s. by pion-quark coupling constants. When we square this 
result and only keep the leading order and the next-to- leading order in l/N^ we finally 
obtain: 

fT^ + Ml = m^9-.T + (-'^^.1. + 2mD.(0)A,-f ) . (4.18) 

For the pion mass we start from Eqs. ( |2.15| ) and ( p.l6|) and expand the inverse pion 
propagator around q'^ = 0: 



2gs n.(0) 



2..(^n.(.) 



+ 0{mt) = . 



(4.19) 



To find in lowest non- vanishing order in mo we have to expand 1 — 2gsIlTT{0) up to 
linear order in mg, while the derivative has to be calculated in the chiral limit, where it 



can be identified with the inverse squared pion-quark coupling constant, Eq. ( 4.15 ). The 
result can be written in the form 



mo 9^ 



qq 



m 2gs 



m 



(4.20) 



Finally one has to expand this equation in powers of l/N^. This amounts to expanding 
9vqq^ which is the only term in Eq. (|4.2CI| ) which is not of a definite order in l/N^. One 
gets: 



m. 



+ 6ml = mn — 



2(0) 



2gs 



mo 



2(0) 
m (jnqq 

2gs 



Jirqq Jnqq 



m 



(4.21) 



It can be seen immediately that the leading-order term is exactly equal to —mo{ipip)^^'^ / fw^^\ 
as required by the GOR relation. Moreover, combining Eqs. ( |2.22|) , (|4.18| ) and (|4.21|) one 
finds that the GOR relation in next-to- leading order, Eq. ( [4.14| ), holds in this scheme. 

However, one should emphasize that this result is obtained by a strict 1/iVc-expansion 
of the various properties which enter into the GOR relation and of the GOR relation 
itself. If one takes /tt and m^r as they result from Eqs. ( [4.17|) and ([4.20|) and inserts them 
into the l.h.s. of Eq. ([4.13|) one will in general find deviations from the r.h.s. which are 
due to higher-order terms in l/Nc- In this sense one can take the violation of the GOR 
relation as a measure for the importance of these higher-order terms [191 . 



4.2 Local selfconsistent scheme 

The proof of the Goldstone theorem in the LSS is very similar to that in the l/N^.- 
expansion scheme. Therefore, we can be brief, concentrating on the steps which are 
different. 

Again we have to show the validity of Eq. (|4.1|) . In the LSS the function 11^^ is given 
by Eq. ( p.6| ), i.e. it differs from to the corresponding function in the l/A'^c-expansion 
scheme (Eq. ( p.l4| )) by the fact, that diagram ^Ili'^'' is (formally) missing. (As we already 
discussed it is implicitly contained in the RPA diagram.) The other diagrams have the 
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same structure as before and we can largely use the results of the previous subsection. 
However, we should keep in mind, that the constituent quark mass is now given by the 
extended gap equation, Eq. ( |3.3| ). Therefore, the r.h.s. of Eq. (O) is different from unity 



in the chiral limit and RPA pions are not massless. This has important consequences for 
the practical calculations within this scheme, which will be discussed in greater detail in 
section |57^ . 

Using Eqs. (^4.8|) , (|4.9| ) and (|4.11| ) as well as Eq. (|3.4| ) we get for the correction terms 
to the pion polarization function 

Y: SUi:\0) = -- = (4.22) 

k=a,b,c 



Hence, together with the modified gap equation (|3.3| ) we find 

2^,11^(0) = 1 - ^ (4.23) 

m 

in agreement with Eq. ( |4.1| ). 

The discussion concerning the regularization procedure can be repeated here. The 
structure of the proof again leads to the conclusion that we have to regularize the quark 
loops in the same way, whereas we have the freedom to choose the regularization for the 
meson loops independently. 

Another important observation is that we, in both schemes, do not need the explicit 
form of the RPA propagators. Da{p) and D^^^p) only need to fulfill Eq. ( [4.7| ). Thus, 
approximations to the RPA propagators can be made as long as Eq. (|4.7| ) remains valid. 

For a non-vanishing current quark mass the pion mass is given by the GOR relation 



(Eq. ([4.13|) ). To linear order in mo this relation holds exactly in the LSS, if we choose 
the appropriate definition of the quark condensate. This will be demonstrated in the 
following. 

For the pion decay constant we follow the same steps as in the 1/iVc-expansion 
scheme to arrive at the following expression: 



_ n^(g)-n^(o) 



Here the modified pion-quark coupling constant is defined as 



. (4.24) 



9ngg = —j^U^=^l ■ (4.25) 



In the chiral limit, 0, the difference ratio on the r.h.s. of Eq. ( [4.24| ) can be replaced 

by the pion-quark coupling constant (Eq. ( [4.25|) . This leads to the Goldberger-Treiman 
relation 

firQ-Kqq = ^ . (4.26) 

Following the analogous steps which led us to Eq. ( [4.20|) we find for the pion mass 



2 

= — ^ + Oim,) . (4.27) 
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Multiplying this with as given by Eq. ( ^.26| ) we get to linear order in mo: 



ml fl = mo ^ . (4.28) 

Obviously this is consistent with the GOR relation (Eq. ( [4.13| )) if the quark condensate is 
given by Eq. (|3.8|) , but not if is given by Eq. ( p. 71) . In Sec. we have seen that within the 
effective action formalism the quark condensate is given by Eq. (|3.8|) . Therefore at this 
point we clearly see that the interpretation of m as a constituent quark mass, which would 
mean that we have to calculate the quark condensate according to Eq. (p.7|), leads to a 
contradiction with the GOR relation. Therefore, in the numerical part, we will calculate 
the quark condensate according to Eq. (|3.8|). 



5 Numerical results at zero temperature 

In this section we present our numerical results at zero temperature. We begin with a 
brief description of the regularization scheme and then discuss peculiarities related to the 
solution of the gap equation in the LSS. After that we study the influence of mesonic 
fluctuations on quantities in the pion sector, thereby focusing on possible instabilities. 
Finally we perform a refit of these quantities within the 1/iVc-expansion scheme and the 
LSS and apply the model to observables in the p-meson sector. 



5.1 Regularization 

Before we begin with the explicit calculation we have to fix our regularization scheme. 
As discussed in Sec. ||, all quark loops, i.e. the RPA polarization diagrams, the quark 
triangles and the quark box diagrams must be regularized in the same way in order to 
preserve chiral symmetry. We use Pauli-Villars-regularization with two regulators, i.e. we 
replace 



(2vr) 



/d k 
^^2^9/(^;^i) ' (5-1) 



with 



/i^' = m^+jA^; Co = l, Ci = -2, Ca = 1 . (5.2) 

Here Ag is a cutoff parameter. 

The regularization of the meson loop (integration over d'^p in Eq. ( p.23| )) is not con- 
strained by chiral symmetry and independent from the quark loop regularization. For 
practical reasons we choose a three-dimensional cutoff A^ in momentum space. In order 
to obtain a well-defined result we work in the rest frame of the "improved" meson. The 



same regularization scheme was already used in Refs. [19, 20 



5.2 Solution of the gap equation in the LSS 

In contrast to the 1 /Ai'c-expansion scheme, where all diagrams are constructed from 



"Hartree" quarks, the LSS is based on the extended gap equation, Eq. (|3l3|). In Sec. p]2 
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this equation was the starting point to find a consistent set of diagrams for the description 
of mesons. In fact, in Sec. |4.2| we have shown, that various symmetry relations, namely 
the Goldstone theorem, the Goldberger Treiman relation and the GOR relation hold in 
this scheme. It is not surprising, that the structure of the extended gap equation was 
needed to prove these relations. So far, all this has been done on a rather formal level. 
This section is now devoted to the explicit solution of the modified gap equation in the 
LSS. We will see, that this cannot be done in a straightforward manner and we are forced 
to a slight modification of the scheme. 

In addition to the Hartree term T^h, Eq. ( p.3|) contains the term which is a quark 
loop, dressed by RPA mesons (see Fig. H). As already pointed out, these RPA mesons 
consist of quarks with the self consistent mass m, which is in general different from the 
"Hartree" mass ttih- Hence, the masses of these mesons are also different from the meson 
masses in the Hartree + RPA scheme. On the l.h.s. of Fig. ^ we have plotted the squared 
masses m^^"^ of the RPA pion (solid) and the RPA cr-meson (dashed) as functions of 
a trial constituent quark mass m'. An important observation is that the pion becomes 
tachyonic, i.e. m^^'^ becomes negative, for quark masses smaller than the Hartree quark 
mass. Strictly speaking this is only the case in the chiral limit, whereas for nonvanishing 
current quark masses, mi°'*^ becomes negative slightly below the Hartree quark mass. A 
similar observation can be made for m^a^"^ , but only for m' much smaller than the Hartree 
mass. This observation of tachyonic RPA mesons is related to the point discussed in 
Sec. |3]^, that the meson-loop term in the effective action (second term of Eq. ( |3.20| )) is 
no longer positive definite. 

Tachyonic RPA mesons lead to a complex correction to the quark self-energy. There- 
fore the solution of the extended gap equation can only be real if it is larger than the 
Hartree mass. Otherwise it must be complex. To investigate this point we plot the dif- 
ference between the l.h.s. and the r.h.s. of Eq. ( |3.3| ) as a function of the (real) trial 
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Figure 9: (Left) Squared pole masses of the pion (solid) and the a -meson (dashed) in 
RPA as functions of a trial constituent quark mass m' in units of the Hartree quark mass. 
(Right) Difference m' — ttiq — S(m') between the l.h.s. and the r.h.s. of the LSS gap 
equation, Eq. ( ^~^ ), as a function of the trial constituent quark mass m'. The real part is 
denoted by the solid line, the imaginary part by the dashed line. 
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quark mass m' . This is shown in the right panel of Fig. ^ The sohd hne denotes the 
real part, the dashed line the imaginary part of m' — mo — S(m'). Obviously, below the 
Hartree quark mass, the self-energy indeed gets complex. Moreover, we see that there is 
no solution of the gap equation for real constituent quark masses. Hence, in principle, 
one should search for solutions of the gap equation in the complex plane. However, this 
would mean that the RPA mesons would also consist of quarks with complex masses. In 
this case, e.g. a reasonable description of the p-meson would be completely impossible, 
because its properties are mainly determined by intermediate pions. 

Therefore, we prefer to perform an approximation, which was introduced in Ref. |[15 |. 
As discussed in Sec. |]2] the symmetry properties of the LSS are not affected by approx- 
imations to the RPA meson propagators which preserve the validity of Eq. ( [4.7|) . The 
authors of Ref. |T^ simply replace the RPA pion propagator in the extended gap equation 



by 



(27r)4 P - m2 + ie 



+ 2igs{2N,Nf)p^I{p) (5.3) 



rriQ 



m 



+ 2igs2N,Nf,p'l{p) 



(5.4) 



and analogously for the a-propagator. The same replacements are performed for the RPA 
meson propagators in the correction terms 5H^^^ to the mesonic polarization diagrams. 
The RPA contribution Hm itself, however, is not changed. In this way the solution of 
the gap equation and the masses of the intermediate mesons remain real. Moreover, in 
the chiral limit the intermediate pions are massless, as one can immediately see from 
Eq. (pD. 

The above replacements would be exact in the Hartree approximation (cf. Eq. ( p.l|) ). 
The authors of Ref. |T5| argue that the correction terms are suppressed because they are 
of higher orders of 1/Nc (beyond next-to-leading order). In the LSS, this is a questionable 
argument because the selfconsistent solution of the gap equation mixes all orders oil/Nc 
anyway. Nevertheless this approximation preserves the validity of the various symmetry 
relations we have checked in Sec. |4.2| . 

In the following we will call this scheme, including the above replacements, the "local 
selfconsistent scheme" although it is strictly speaking only an approximation to the LSS 
as it was originally introduced in Sec. Ol 



5.3 Meson-loop effects on quantities in the pion sector 

In this subsection we want to study the influence of mesonic fluctuations on the quark 
condensate, the pion mass and the pion decay constant, both within the l/A^^c-expansion 
scheme and within the LSS. Since the strength of the fluctuations is controlled by the 
meson cutoff A^, we first keep all other parameters fixed and investigate how the above 
quantities change, when Km is varied. For the l/A'^c-expansion scheme this has been done 
in more detail in Ref. [^. Later, in the next subsection, we will perform a refit of the 
parameters to reproduce the empirical values of {'ipip), and 

Our starting point is the Hartree -|- RPA scheme, which corresponds to Am = 0. 
Here we obtain a reasonable fit {{ijjip)^^^ = -2 (241.1 MeV)^ mi°^ = 140.0 MeV and 
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Figure 10: The ratios ml/ ml^'^'^ (solid), fH ff"^^ (dashed), (dashed- dotted) 

and the combination —'mo{'ip'ip)/m'lf^ (dotted) as a function of the meson loop cutoff Am- 
Left: 1/Nc- expansion scheme. Right: LSS 



/^°) = 93.6 MeV) with the parameters = 800 MeV, ^.A^ = 2.90 and mo = 6.13 MeV. 
These parameters correspond to a relatively small "Hartree" constituent quark mass of 
260 MeV. 

Now we turn on the mesonic fluctuations by taking a non-zero meson cutoff Am- All 
other parameters are kept constant at the values given above. The resulting behavior of 

is displayed in Fig. 10. The left 



ml, f^ and the quark condensate as a function of Am 
panel corresponds to the l/A^^c-expansion scheme, the right panel to the LSS. As one can 
see, in both schemes the mesonic fluctuations lead to a reduction of f-,, (dashed lines) 
while m,r (solid) is increased. At smaller values of Am the absolute value of the quark 
condensate decreases but goes up again for Am ^ 900 MeV. This is also an effect which 
is found in both schemes. 

In the Hartree + RPA scheme the quantities ml^^\ f^^^^ and {■iljip)^^\ are in almost 



perfect agreement with the GOR relation, Eq. ( [4.13|) . As discussed in Sec. ^3], the 1/Nc- 
expansion scheme is consistent with the GOR relation up to next-to-leading order in 
1/Nc, but the relation is violated by higher-order terms. We therefore expect a less 
perfect agreement in this scheme, becoming worse with increasing values of Am- In the 
LSS, on the other hand, the quantities should be in good agreement with the GOR (see 
Sec. 13). 

These expectations are more or less confirmed by the results. In Fig. |T0|, the ratio 
of the r.h.s. and the l.h.s. of Eq. ( [4. 13] ) is displayed by the dotted lines. In the 1/Nc- 
expansion scheme (left panel) the relation holds within 30% for Am < 900 MeV. However, 
when the meson cutoff is further increased the deviation grows rapidly. This indicates that 
higher-order corrections in 1/Nc become important in this regime and this perturbative 
scheme should no longer be trusted. In the LSS the agreement with the GOR is almost 
perfect. 

In Fig. |10| the various curves are only shown up to Am = 1250 MeV for the 1/Nc- 
expansion scheme and Am = 950 MeV for the LSS. For larger values of Am a second, 
unphysical, pole with a residue of the "wrong" sign, emerges in the pion propagator. 
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This would correspond to an imaginary pion-quark coupling constant and hence an imag- 
inary pion decay constant. Upon further increasing Km the two poles merge and finally 
disappear from the positive real axis. 

For the l/A^^c-expansion scheme this has been discussed in more details in Ref. In 
that reference we suggested that the instabilities of the pion propagator might indicate 
an instability of the underlying ground state against mesonic fiuctuations. In fact, it has 
been claimed by Kleinert and Van den Bossche |^ that there is no spontaneous chiral 
symmetry breaking in the NJL model as a consequence of strong mesonic fiuctuations. 
Although this cannot be true in general if the strength of the mesonic fiuctuations is 



controlled by an independent cutoff parameter Am [|T9[5 this phenomenon might very 
well occur for large values of A^- In other words: There could be some kind of "chiral 
symmetry restoration" at a certain value of the parameter A^- 

Clearly, this could not be studied within the l/A'^c-expansion scheme where the mesonic 
fluctuations are built perturbatively on the Hartree ground state. In the LSS, however, 
where we encounter the same type of instabilities in the pion propagator, this question 
can be investigated more closely. To that end we consider the effective action Eq. ( |3.20| ), 
which describes the energy density of the system. It is explicitly given by 



Tim') = -AiN.Nf f /4ln( '", ^ J + 



^9s 



-\ j - '^9s^M) + 31n(l - 2gsYlM)} + const. . (5.5) 

The irrelevant constant can be chosen in such a way that r(0) = 0. The positions of the 
extrema of V{m') correspond to the solutions of the gap equation ( |3.3|) . In particular, 
the vacuum expectation value m is given by the value of m' at the absolute minimum of 
r. Note that, according to Eq. (|3.8| ), m is proportional to the quark condensate, i.e. to 
the order parameter of chiral symmetry breaking. Hence, for a given value of Km-, chiral 
symmetry is spontaneously broken, if the absolute minimum of F is located at a non-zero 
value of m' and it is unbroken ( "restored" ) otherwise. 

We perform the calculations in the chiral limit .[| Our results for F(m') as a function of 
m! /mu for different values of A^ are displayed in Fig. |TT|. For Km = we flnd of course 
the minimum at m' = itlh = 260 MeV, while there is a maximum at m' = 0. If there was 
indeed a "phase transition" due to mesonic fluctuations, this maximum should eventually 
convert to a minimum when Am is increased. In fact, for Am ^ 900 MeV the results 
seem to point in this direction: In this regime the constituent quark mass m is reduced to 
about 30% of the Hartree mass. At the same time the "bag constant" B = F(0) — F(m) 
decreases from 48.7 MeV/fm^ at Am = to 0.8 MeV/fm^ at Am = 900 MeV. However, 
upon further increasing Am, both m and B go up again. In particular, the point F(0) 
always remains a local maximum: In the LSS we do not observe a "phase transition" due 
to strong mesonic fluctuations. 

Here we should remark, that, according to the conjecture by Kleinert and Van den 



Bossche 22 , the mesonic fluctuations do not restore the trivial vacuum in the NJL model, 



TTo be precise, we proceed as follows: Starting from the parameters given above, we keep the Hartree 
constituent quark mass, itih ~ 260 MeV, fixed, while mg is reduced from 6.1 MeV to zero. Therefore the 
coupling constant is slightly enhanced from .g^A^ = 2.90 to (?sA^ — 2.96. 
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Figure 11: Effective action r(m') as a function ofm'/mH in the LSS for different values 
of the meson cutoff Am: MeV (solid), 300 MeV (long- dashed), 500 MeV (dotted), 
900 MeV (dashed- dotted) and 1200 MeV (short- dashed). 



but lead to a so-called pseudo-gap phase. (See also Ref. for a critical discussion of 
this article.) In that phase the quarks still have a non-vanishing constituent mass, if the 
latter is identified with the vacuum expectation value of the modulus of the scalar field 
$ (cf. Sec. |3.3|) . Nevertheless chiral symmetry is not broken as the phase of the $ field 
is strongly fluctuating. An analogous phenomenon is well known from strong-coupling 
superconductors above Tc [|^, 0, where Cooper pairs are formed, but do not condense. 
Obviously our above investigations, which focused on a the change of m assuming a 
uniform phase factor, cannot exclude a transition into a phase of this type. Here more 
refined investigations are needed to give a conclusive answer. 

Another type of vacuum instability which is caused by unphysical poles of the RPA 
meson propagators has recently been discussed by Ripka Here "unphysical" means 
that these poles are located in regions of the complex plane, where they are forbidden by 
microcausality. Ripka stated, that they are induced by the regulator scheme, in his case a 
4-momentum cutoff or a Gaussian form factor. In fact, the RPA meson propagators have 
this unphysical feature for most of the known regulator schemes, such as proper-time reg- 
ularization, subtracted dispersion relations, dimensional regularization or, as mentioned 
above, a 4-momentum cutoff. A 3-momentum cutoff and Pauli-Villars-regularization in 
the form we use it (cf. appendix P) are exceptions. On the other hand, due to Pauli- 
Villars regulators, the imaginary part of the quark loops can have the wrong sign in some 
kinematical regions and we cannot rule out that the instabilities we find for the pion 
propagator are related to this. This supposition is corroborated by the fact that these 
instabilities could be traced back to the imaginary part of the diagram (see Fig. 
which has the "wrong" sign and which becomes large at large values of Am [|1^]- Further 
investigations are needed, however, to clarify this point. 

Recently a second (unphysical) pole in the pion propagator has also been found in 
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a non-local generalization of the NJL model The calculations indicate that these 

instabilities could probably be removed by including vector and axial vector intermediate 
states. This point is certainly worth a closer examination. In any case, at least in the 
1/iVc-expansion scheme we found [^] that with a reasonable fit of all parameters we are 
far away from the region where these instabilities occur. We will come back to this point 



in Sec. 5.5 



5.4 Parameter fit in the pion sector 

In the previous subsection we did not change the parameters which were determined in 
the Hartree + RPA scheme by fitting fi^\ mi°^ and {ipip)'^'^^ Of course, if one wants 
to apply the model to describe physical processes a refit of these observables should be 
performed including the mesonic fluctuations. In Ref. this was already done for 



the 1/iVc-expansion scheme and we will now try to perform an analogous fit within the 
LSS. Of course, by fitting the above three observables, we cannot conclusively fix the 
five parameters of our model, gs-, Qv, Ag, Am and rriQ. Therefore we try to proceed in a 
similar way as in Ref. pO[: For various values of Am we fix the scalar coupling constant 



gs, the current quark mass mo and the quark- loop cutoff Ag to reproduce the empirical 
values of the pion mass, /^r and (ipip). Then, in the next subsection, we will try to 
fix the two remaining parameters, i.e. the vector coupling constant Qy and the meson 
cutoff Am, by fitting the pion electromagnetic form factor in the time-like region, which 
is related, via vector meson dominance, to the p-meson propagator. Roughly speaking, 
this amounts to fitting the p-meson mass and its width. Since in our model the latter is 
due to intermediate RPA pions, we decided to fix the empirical value of not m^, in 



order to get the correct threshold behavior. In Ref. ||2^ we found for the l/A'^c-expansion 
scheme that the deviation is about 10%. As we will see below, in the LSS the difference 
is somewhat larger. 

Of course, the p-meson can only be described reasonably if the unphysical gg-threshold 
lies well above the peak in the p-meson spectral function, i.e. the constituent quark mass 
m should be larger than about 400 MeV. For that reason we try to increase the constituent 
quark mass as much as possible. Here we have some freedom as the empirical value of the 
quark condensate is not known very precisely. (Its absolute value is probably less than 
2(260 MeV)^, which corresponds roughly to the upper limit extracted in Ref. ||3^ from 
sum rules at a renormalization scale of 1 GeV. Recent lattice results give (ipip) = -2((231± 
4 ± 8 ± 6) MeV)^ 0-) On the other hand, since the correction term in the LSS gap 
equation, Eq. ( p.3|) , contributes negatively to m, it is much more difficult in the LSS to 
obtain sufficiently large quark masses than in the l/A^^c-expansion scheme. 

Our results for the LSS are given in Table 0. For comparison we also summarize the 
results obtained in Ref. |^ within the 1/A^c-expansion scheme (Table |1|). In both tables 
we list five parameter sets (corresponding to five different meson cutoffs Am), together 
with the constituent quark mass m, the values of m^r, /jr and (ipip) and the corresponding 
RPA quantities. In the LSS the "RPA quantities" are calculated with the constituent 
quark mass m in order to represent the properties of the intermediate pion states. For 
completeness we also give the value of the Hartree mass mn in Table ^ and the value of the 
quark condensate according to Eq. (^.7|). We also show the ratio —mo{ijjip)/iTi^f^, which 
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Am I MeV 


0. 


300. 


500. 


600. 


700. 


Ag / MeV 


800. 


800. 


800. 


820. 


852. 


mo / MeV 


6.13 


6.40 


6.77 


6.70 


6.54 


9sA^ 


o on 


o.U / 


O Ad 

o.4y 


o. / U 


A A p. 

4. io 


m 1 MeV 


ZOU. 


oU4. 


oyo. 


44d. 


boU. 


m^r / MeV 


140.0 


140.0 


140.0 


140.0 


140.0 


/ MeV 


140.0 


143.8 


149.6 


153.2 


158.1 


jr 1 MeV 


93.6 


100.6 


111.1 


117.0 


126.0 


1 MeV 


93.6 


93.1 


93.0 


93.1 


93.4 


/ MeV^ 


-2(241.1)3 


-2(249.3)3 


-2(261.2)3 


-2(271.3)3 


-2(287.2)3 


{^f) 1 MeV^ 


-2(241. If 


-2(241.7)3 


-2(244.1)3 


-2(249.5)3 


-2(261.4)3 




1.001 


1.007 


1.018 


1.023 


1.072 



Table 1: The model parameters (Am, Ag, mo and gs) and the resulting values ofrriT,, f.^ 
and {ijjip) (together with the corresponding leading-order quantities), the constituent quark 
mass m in the l/N^- expansion scheme. The ratio —mo{il!tp)/m'^f^, is also given. 



Am I MeV 


0. 


300. 


500. 


600. 


700. 


Ag / MeV 


800. 


800. 


810. 


820. 


835. 


mo / MeV 


6.13 


6.47 


7.02 


7.30 


7.90 


9sAl 


2.90 


3.08 


3.44 


3.71 


4.52 


ran / MeV 


260. 


305. 


390. 


450. 


600. 


m / MeV 


260. 


278.2 


320.0 


355.7 


468.4 


ml°^ / MeV 


140.0 


139.9 


140.0 


139.7 


140.0 


m^ / MeV 


140.0 


145.1 


156.3 


164.5 


182.7 


fi°^ / MeV 


93.6 


96.7 


103.6 


108.4 


120.0 


U / MeV 


93.6 


93.2 


92.9 


92.9 


92.8 


(il^i))' / MeV3 


-2(241.1)3 


-2(244.7)3 


-2(254.3)3 


-2(261.9)3 


-2(277.3)3 


(iljip) / MeV3 


-2(241.1)3 


-2(241.7)3 


-2(246.2)3 


- 2(250.8)3 


-2(260.9)3 




1.001 


1.001 


1.006 


1.01 


1.02 



Table 2: The same as in Table |7| for the LSS. The quantity (ipip)' denotes the quark 
condensate calculated according to Eq. ^3. 

would be equal to 1 if the GOR relation was exactly fulfilled. Note that the deviations 
in the 1 /A^c-expansion scheme are less than 10% (for A^/ < 600 MeV even less than 3%), 
indicating that higher-order corrections in 1/A"c are small. In the LSS the deviations are 
considerably smaller, as already discussed in the previous subsection. 

In both schemes we find that the constituent quark mass increases with an increasing 
meson cutoff Am- In the l/A'c-expansion scheme for A^ > 500 MeV the value of m is large 
enough to shift the gg-threshold above the p-meson peak. Besides, it turns out that we can 
only stay below the limit of — 2(260Mel^)3 for the quark condensate and simultaneously 
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reproduce the empirical value of if the cutoff is not too large (Am ^ 700 MeV). In 
the LSS the region of values for where on one hand the constituent quark mass is 
large enough and on the other hand the quark condensate stays below the limit is much 
more narrow. This can be seen from the values listed in Table ^ For a meson cutoff of 
Am = 600 MeV m is still too small and for Am = 700 MeV the quark condensate lies 
already slightly above the limit. The reason for this is obvious: In the LSS m and the 
quark condensate are directly related by Eq. (|3.8|) and therefore the mesonic fluctuations 
which lower the quark condensate also decrease the constituent quark mass. In the l/N^- 
expansion scheme, on the contrary, the meson loop effects only contribute to the quark 
condensate and lower its value whereas m is kept fixed at its Hartree value. 



5.5 Description of the /o-meson 

As already pointed out, the parameter fit in the pion sector was not complete. It is 
clear, e.g., that the meson- loop cutoff Am cannot be determined just by fitting the pion 
mass, and (ipip), since these observables can already be reproduced in the Hartree 
+ RPA scheme, i.e. without any meson-loop effects. We only found an upper limit of 
Am ^ 700 MeV in both schemes (see Tables |l| and H). In Ref. |^ we therefore fixed the 



remaining parameters gv and Am for the l/A^^c-expansion scheme in the p-meson sector. In 
this subsection we want to give a short summary of these results and then try to perform 
a similar fit for the LSS. 

According to Eqs. ( p.l4|) and (|3.6| ), the polarization function of the p-meson reads: 



n^'^'"^(g) = n^'^'"^(g) + J2 sul^^ """'"^q) ■ (5.6) 

k 

Here k runs over {a,b,c,d} in the l/Ai"c-expansion scheme and only over {a,b,c} in the 
LSS. Because of vector current conservation, the polarization function has to be transverse, 

g^Hf '«"(?) = q.fl';^'^\q) = 0. (5.7) 

With the help of Ward identities it can be shown that these relations hold in both schemes, 
if we assume that the regularization preserves this property. This is the case for the 
Pauli-Villars regularization scheme, which was employed to regularize the RPA part Hp. 
Together with Lorentz covariance this leads to Eq. { ^■S\) for the tensor structure of Hp. 
On the other hand, since we use a three-dimensional sharp cutoff for the regularization of 
the meson loops, the correction terms ^Hp'^^ are in general not transverse. However, 
as mentioned in Sec. |5.1|, we work in the rest frame of the p-meson, i.e. g* = 0. In this 



particular case Eq. (|5.7|) is not affected by the cutoff and the entire function Hp can be 



written in the form of Eq. 

Il^'^'^'ig) = %{q)T'^''6,, = (Hp(g) + J] <5HW(g)) T^'^ 5,, , (5.8) 

k 

i.e. instead of evaluating all tensor components separately we only need to calculate the 
scalar functions Hp = — l/Sgfpi.H'^'' and SU^^^ = —l/Sg^uSH^p^^'^. 
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Figure 12: Contributions to the pion electromagnetic form factor in the 1/ N^.- expansion 
scheme. The propagator denoted by the curly line corresponds to the 1/ N^.- corrected rho- 
meson, while the double lines indicate RPA pions and sigmas. 



A second consequence of vector current conservation is, that the polarization function 
should vanish for = 0. For the correction terms this is violated by the sharp cutoff. 
We cure this problem by performing a subtraction: 



(5.9) 



Note, however, that already at the RPA level a subtraction is required, although the RPA 
part is regularized by Pauli-Villars. This is due to a rather general problem which is 
discussed in detail in App. ^. 

In Ref. I^Ol we have fixed and Am in the l/A^^c-expansion scheme, by fitting the 
pion electromagnetic form factor, F^^lq), in the time-like region, which is dominated by 
the p-meson. The diagrams we included in that calculations are shown in Fig. The 
two diagrams in the upper part correspond to the standard NJL description of the form 



factor [34] if the full p- meson propagator (curly line) is replaced by the RPA one. Hence, 
the first improvement is the use of the 1/A^c-corrected p- meson propagator in the 1/Nc- 
expansion scheme. Since, in the standard scheme, the photon couples to the p-meson via 
a quark-antiquark polarization loop, in the l/A'c-expansion scheme we should also take 
into account the l/A^^c-corrections to the polarization diagram for consistency. This leads 
to the diagrams in the lower part of Fig. On the other hand the external pions are 
taken to be RPA pions (i.e. mass rrin^ and pion-quark-quark coupling constant gnqg)- 
This is more consistent with the fact that the p-meson is also dressed by RPA pions and, 
as discussed above, we have fitted m^^ to the experimental value. 

The numerical results for l-F^rp as a function of the center-of-mass energy squared are 
displayed in the left panel of Fig. |T3|, together with the experimental data |^ . The the- 
oretical curve was calculated with a meson cutoff of Am = 600 MeV, a vector coupling 
constant g^ = l.Ggs and the other parameters, Ag, gs and mo as listed in Table |T]. This 



roughly corresponds to a best fit to the data [Q. Since we assumed exact isospin sym- 
metry we can, of course, not reproduce the detailed structure of the form factor around 
0.61 GeV^, which is due to p-tu-mixing. The high-energy part above the peak is somewhat 
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Figure 13: The pion electromagnetic form factor (left panel) and the irir-phase shifts in 
the vector-isovector channel (right panel) for Am = 600 MeV and = l.Qgs- The other 
parameter values are taken from Table The data points are taken from refs. ^3d{J and 
respectively. 



underestimated, mainly due to the sub-threshold attraction in the p-mesons channel be- 
low the gg-threshold at s = 0.80 GeV^. Probably the fit can be somewhat improved if we 
take a slightly larger meson cutoff, but we are not interested in fine-tuning here. 
A closely related quantity is the charge radius of the pion, which is defined as 

dq" 

With the above parameter set we obtain a value of (r^)^/^ = .61 fm. It lies slightly below 
the experimental value, (r^)^^^ = (0.663 ± 0.006) fm . 



(-^) = 6^. • (5-10) 
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Figure 14: Diagrams contributing to the tttt -scattering amplitude: Quark box diagram 
(left) and s-channel p-meson exchange (right). 

One can also look at the vrvr-phase shifts in the vector-isovector channel. We include 
the diagrams shown in Fig. [1^, i.e. the s-channel p-meson exchange and the direct tttt- 
scattering via a quark box diagram. The latter has to be projected onto spin and isospin 1, 
which is a standard procedure. (For example, the analogous projection onto spin and 
isospin can be found in Refs. ^9|, ^.) The result, together with the empirical data 
3^, is displayed in the right panel of Fig. |TB|. Since the main contribution comes from 



the s-channel p-meson exchange, it more or less confirms our findings for the form factor: 
below the p-meson peak a good of the data is obtained while, at higher energies, where 
gg-threshold effects start to play a role, we slightly overestimate the data. 
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Figure 15: Real part of the p-meson polarization function Yip multiplied by 2gy = 
17.6 GeV'"^ as a function of the energy -s/s in the rest frame of the meson. The dashed 
line corresponds to the 1/Nc- expansion scheme with Am = 600 MeV, the solid line to 
the LSS with Am = 700 MeV. The other parameters are given in Table |I| and Table 
respectively. 



Let us now turn to the LSS. As already discussed in the last paragraph of Sec. 5.4 



there is not much room to vary the meson cutoff Am in this scheme: On the one hand Am 
is restricted to values < 700 MeV by the fit to and (ipip) (see Table |]). On the other 
hand we only have a chance to get a realistic description of the p-meson if the constituent 
quark mass m is larger than at least 400 MeV. To achieve this, the meson cutoff cannot 
be much smaller than 700 MeV. This means, Am is more or less fixed to this value and 
the only remaining parameter is the vector coupling constant g^. 

It turns out, however, that with Am 700 MeV we run already into instabilities in the 
p-meson channel. These instabilities are of the same type as the instabilities in the pion 
channel, (see Sec. ^.3|) , but unfortunately emerge already at lower values of A^. This 
can be seen in Fig. |T3| where the real part of the of the p-meson polarization function Up 
multiplied by 2gv is plotted as a function of the energy y/s in the rest frame of the meson. 
The LSS result corresponds to the solid line. For comparison we also show this function 
in the 1 /Ac-expansion scheme, using the 'best-fit parameters' given above (dashed curve). 

According to Eqs. ( |2.15| ) and ^TS| the function 2g^ ReUp has to become equal to 1 for 
^/s ~ rup, crossing the line 2gy ReUp = 1 from below. This is obviously the case in the 
1 /A"c-expansion scheme. In this scheme, i^ellp is a rising function and the above condition 
can be easily fulfilled with the appropriate choice of g^. The situation is quite different in 
the LSS. Here Rellp has a maximum at ^/s ~ 740 MeV and then steeply drops. Hence, 
if gy is too small, the equation 2gy Rellp = 1 has no solution at all (see Fig. |15|). On 
the other hand, for large values of g^ we get a "phsyical" solution at lower energies and 
an "unphysical" solution at higher energies. It is clear that none of these two scenarios 
would lead to a realistic description of the p-meson. 
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Figure 16: RPA contribution 2gyllp (solid) and the various correction terms to 2g^I\.p: 
2g^5Ii^f' (dashed), 2g^{5Ii^p'^ + dli^f^) (dotted) and2gy6Yi^f' (dashed- dotted). For all con- 



tributions we performed a subtraction, such that they vanish at a/s = 0. The left panel 
corresponds to the 1/Nc- expansion scheme, the right panel to the LSS. The model param- 
eters are the same as in Fig. [7^. 



One might wonder, why the results in the 1 /A'^c-expansion scheme and in the LSS are 
so different. To answer this question we have separately plotted the various contributions 
to the polarization function in Fig. The left panel corresponds to the results in 
the 1/A^c-expansion scheme, the right panel to the LSS. One immediately sees that the 
unphysical behavior in the LSS is due the sum of the diagrams (5n?^ and 6Ui'^ (dotted), 
which is the only negative contribution. In the 1/iVc-expansion scheme these diagrams 
behave very similar. However, in this scheme their contribution is almost cancelled by 
the contribution of diagram ^LIp"''*, which is not present in the LSS. 

We should note that diagram ^LIp"'', which describes the two-meson intermediate state 
is well-behaved in both schemes. On the other hand the momentum dependence of all 
other diagrams is a pure quark effect, which could be related to the imaginary part 
of these diagrams above the (unphysical) two-quark threshold via dispersion relations. 
Hence, if we could manage to further push up the constituent quark mass, the momentum 
dependence of these contributions should become smaller and the instabilities should 
eventually vanish. Perhaps this is possible if further intermediate mesons, like p and Oi 
are included in the model. 



6 Quark condensate at T 7^ 

It is expected, that at sufficiently large temperatures chiral symmetry, which is spon- 
taneously broken in vacuum, gets restored. The quark condensate as an order parame- 
ter of chiral symmetry is well suited to study indications for (partial) chiral symmetry 
restoration. At low temperatures model independent results for the changes of the quark 
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condensate can be obtained from considering a gas of pions, which are the hghtest par- 
ticles and therefore the main degrees of freedom in this range. Approaching the phase 
transition we have to rely on model calculations or lattice data because we do not have 
any fundamental knowledge of the quark condensate at higher temperatures. Most of the 
results show a phase transition at a temperature of Tc ~ 150 MeV. 

Among others, the NJL model has been used to examine the behavior of the quark 
condensate as a function of temperature. Most of these investigations have been performed 
in the mean- field approximation |]^, ^, ^ 0. There one finds a second order phase 
transition with Tc ~ 150 — 200 MeV. However one has to mention that these calculations 
suffer from the severe problem that the thermodynamics is generated exclusively by a 
gas of quarks. One consequence is that the low temperature behavior, which is driven 
mainly by pions, is completely missed. Although we cannot by-pass the fundamental 
problem of lack of confinement in the NJL model which in any case leads to the existence 
of a quark gas at non-zero temperature, we can hope to improve the situation at least 
at low temperatures via inclusion of mesonic degrees of freedom in a calculation beyond 
mean-field. 

Therefore we begin with a closer look at the low-temperature behavior of the quark 
condensate at T 7^ 0. After that we will discuss our numerical results within the 1/Nc- 
expansion scheme and within the LSS. 



6.1 Low-temperature behavior 

In the chiral limit and at vanishing baryon density a strict low-temperature expansion in 



chiral perturbation theory leads to the following expression for the quark condensate ^] 



(^^)t = (^^)(i-^-^ + ...) • (6.1) 

Here (ipip) denotes the quark condensate at zero temperature. The T^-term represents the 
contributions from a pure pion gas, whereas the higher-order terms are due to interactions 
between the pions. It has been shown that the T^- and the T^-term of this expansion 
are model independent results which follow from chiral symmetry alone. Thus in principle 
every chirally symmetric model, including the NJL model, should reproduce these terms. 
However, as is of the order \/W^, we can see that they are of the order 1/A^c and 
respectively. So a mean-field calculation, which corresponds to a restriction to leading in 
1/Nc, will not be able to reproduce these terms |23]. Indeed, NJL model calculations in 



mean-field show a much more fiat behavior at low temperatures |T^, : 



mP = { 1 - ^^?^e^f + ...). (6.2) 



TT 



3/2 (^^)(0) 



Extending the calculations to next-to- leading order in l/N^. will allow us to reproduce the 
T^-term. This will be demonstrated in the following. 

Our calculations at non-zero temperature are performed within imaginary time for- 
malism. Basically this amounts to replacing the energy integration in the various n-point 
functions by a sum over Matsubara frequencies. The explicit expressions are listed in 



37 



App. ^ As there exists a preferred frame of reference in the heat bath, all dynamical 
quantities depend separately on energy and three-momentum. Hence in the following, a 
finite-temperature RPA propagator, for instance, will be denoted as DMifjJ.v)- Foi^ scalar 
quantities, like masses or condensates at non-zero temperature we use a suffix T in order 
to distinguish them from the analogous quantities in vacuum (cf. Eqs. and (|6.2| )). 

In analogy to the vacuum expressions (Eqs. ( p.5|) and ( p.22| )) the quark condensate in 
next-to- leading order of the l/A'^c-expansion scheme is given by 

/7/\ /7 /\(o) , x/7 /\ rriT-mo D^{0,0)At , . 

(Wt = W)t +'^(Wt = • (6.3) 

As shown in Eq. (|6.2|) , the leading-order term (ipip)^^ does not contribute to the change 
of the quark condensate to order T^. Similarly, thermal effects in the a- meson propagator 
can be neglected at low temperatures. Therefore we only need to consider the temper- 
ature dependence of A^. If standard techniques are used the sum over the Matsubara 
frequencies in Eq. ( p.5| ) can be converted into a contour integral f^ : 



= AtN,NjmT^J-0^^l^-^^{D^{z,p) {31(0,0) -3{z'- 

+D^iz,p) (2 /(z,p) + /(0,0)-(z2-pl-4m2.) Kiz,p))^ .(6.4) 

At low temperatures, the main contribution to the temperature-dependent part of this 
integral comes from the lowest lying pion pole, as the other contributions are exponentially 
suppressed. In the chiral limit we can therefore approximate this part for low temperatures 
by 

A.-A = 4iV.iV,m/|?3^{^^}. (6,5) 
This integral can be evaluated analytically and we obtain: 

At-A = m- . (6.6) 

The last step is to realize that in the chiral limit the vacuum cr-meson propagator can be 
expressed through the leading-order pion decay constant as 

^^(0) = (6-7) 



see Eqs. ( p.4|) and ( B.9 )). We finally obtain for the quark condensate in next-to-leading 



order at low temperatures: 

(V^V^)t = (^^)-(V^V^)(°)^ (6.8) 

Comparing this with the chiral perturbation theory result, Eq. (|6.1|) , we see that we can 
in principle reproduce the T^-term. Note, however, that the coefficient is given by the 
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quark condensate and the pion decay constant in leading order in l/A'c, according to a 
strict expansion of Eq. (|6.1| ) up to next-to- leading order in l/Nf.. The physical reason for 
this behavior is the fact that the 1/Nc corrections to the quark condensate correspond to 
fluctuating RPA-mesons and hence the thermal corrections at low temperatures are due 
to thermally excited RPA pions in this model. 

For the LSS, a similar result has been derived in Ref. [^. In the chiral limit the 
authors find 



1 



2(0) 



(6.9) 



Here Z^*-"^ is understood as the RPA-pion decay constant, Eq. ( |B.9D , but evaluated at the 
quark mass m, which follows from the LSS gap equation, Eq. ( p.3|) . This corresponds to the 
fact that in the LSS the thermal corrections to the quark condensate at low temperatures 
are due to RPA pions which consist of LSS quarks. 



6.2 Numerical results within the 1/A^c-expansion scheme 

Our numerical results for the temperature behavior of the quark condensate within the 
1 /A^c-expansion scheme are displayed in Fig. The r.h.s. corresponds to a realistic 
parameter set with mi°^ = 140 MeV (Table I with Am = 600 MeV), the l.h.s. to the 
chiral limit. The solid lines indicate the results obtained in next-to- leading order. For 
comparison we also show the leading order (dashed line) and the pure pion gas result 
(dotted). 

We begin our discussion with the chiral limit. At low temperatures (T < 100 MeV) 
our results show the behavior discussed in the previous subsection: The next-to-leading 
order result is in very good agreement with the pion gas result (Eq. |6.8|) , whereas the 
leading-order result remains almost constant. Therfore in this regime the extention of the 
NJL model to next-to-leading order in l/A'c leads to a considerable improvement. Since 
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Figure 17: Quark condensate as a function of temperature, normalized to the vacuum 
value, in the chiral limit (left) and with mir*-* = 140 MeV (right). Leading order in l/N^ 
(dashed line), next-to-leading order (solid line) and free pion gas (dotted line). 
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the unphysical quark degrees of freedom, which are in principle always present in the 
NJL model, are exponentially suppressed, the system is dominated by the (physical) pion 
degrees of freedom, which come about in next-to-leading oder. 

However, because of the much larger degeneracy factor (24 as compared to 3) we 
cannot avoid that effects due to thermally excited quarks become important at some 
temperature. In our present calculation this happens at about T ~ 100 MeV. In a free 
gas approximation, this roughly corresponds to the temperature, at which the quark 
pressure becomes equal to the pion pressure. 

At this point one might raise the question about the physical meaning of quark effects 
at these temperatures. In nature, quark degrees of freedom can only be excited above 
the deconfinement phase transition. In the NJL model there is no confinement and hence 
no deconfinement transition. However, lattice calculations indicate that the decon- 



finement phase transition at finite temperature coincides with the chiral phase transition. 
One should therefore compare the temperature at which thermally excited quarks become 
important with the critical temperature for the chiral phase transition. From the above 
point of view, quark effects below the phase transition are either invisible or unphysical. 
On the other hand, at least close to the phase transition one might relax this strict po- 
sition. In this regime one might think of a resonance gas with many degrees of freedom, 
which could be effectively described by a quark gas ( "quark-hadron duality"). 

Unfortunately, as already pointed out in Sec. ||, the perturbative treatment of the 
mesonic fluctuations does not allow for a description of the chiral phase transition. Al- 
though the quark condensate vanishes at T ~ 200 MeV, this does not correspond to a 
true phase transition. (Note that the slope of the curve does not diverge at this point.) 
In any case, the applicability of the perturbative expansion scheme probably breaks down 
much earlier. Therefore we cannot give a definite answer to the question whether the 
thermally excited quarks become important near the phase transition or much below. 

Our results with mo 7^ are shown on the r.h.s. of Fig. [l^. Since the RPA pions are 
now massive and therefore exponentially suppressed, the quark condensate as a function 
of T stays much more fiat than in the chiral limit. Nevertheless, at low temperatures pions 
can still be most easily excited as they are the lightest particles. Therefore the next-to- 
leading order result (solid line) can be approximated quite well albeit not perfectly by the 
pure pion gas result (dotted) in this regime. The latter was calculated from the pressure 
Pt, of a massive pion gas as 

= m + m''^^^ (6.10) 

which can be easily derived with the help of the GOR relation. 

Quark effects become important at almost the same temperature as in the chiral limit, 
at T ~ 100 MeV. 

6.3 Local selfconsistent scheme 

Let us now compare the results of the previous subsection with the analogous calculations 
in the LSS. A study of the temperature dependence of the quark condensate within the 
LSS can also be found in Ref. |2^. Here we restrict ourselves to the chiral limit. 
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Figure 18: Left: Quark condensate in the chiral limit as a function of temperature, nor- 
malized to the vacuum value, Hartree approximation (dashed), local self consistent scheme 
(solid) and free pion gas (dotted). Right: Thermodynamic potential per volume as a 
function of the constituent quark mass in the self consistent scheme with T = 163.9 MeV 
(dotted), T = 164.5 MeV (dashed) and 165.3 MeV (solid). 



Our results are shown on the Lh.s. of Fig. 18. The calculations have been performed 



using the parameters of Table y for Am = 700 MeV, but mo = 0. As discussed in 
Sec. |0| , at low temperatures the model behaves again like a free pion gas (dotted line). 
Deviations from this behavior become visible at T ~ 100 MeV which is quite similar to 
our observations in the l/A^^c-expansion scheme. 

In contrast to the 1/A^c-expansion scheme, the treatment of the mesonic fluctuations in 
the LSS also allows an examination of the phase transition. With the present parameters it 
takes place at = 164.5 MeV which is considerably lower than in Hartree approximation, 
where we have = 266.1 MeV. Note, however, that about one third of this reduction 
can be attributed to the fact that the constituent quark mass m = 468.4 MeV in the LSS 
is lower than the corresponding Hartree mass mn = 600.0 MeV. For mn = 468 MeV we 
would get a critical temperature of about 236 MeV in Hartree approximation. (It is also 
interesting to note, that the critical temperature in the LSS calculation almost coincides 
with the critical temperature T^^^ = 164.4 MeV one obtains in Hartree approximation 
for the parameters fitted in the RPA, i.e. m = rriH = 260 MeV.) 

Whereas in Hartree approximation the phase transition is of second order, in the 
LSS the system undergoes a first-order phase transition, as already reported in Ref. p3[| . 
This can be inferred from the thermodynamic potential u, which is displayed on the 
r.h.s. of Fig. |18] for different temperatures as a function of the constituent quark mass 
m'. At T = 164.5 MeV one can clearly identify two degenerate minima at m' = and 
m' 7^ 0, corresponding to a first-order phase transition at that temperature. One can ask 
whether this phenomenon depends on the strength of the mesonic fluctuations which can 
be controlled by the cutoff A^. Varying this parameter we find that the discontinuity 
decreases with decreasing Am, but even for very small values of the cutoff we encounter 
a first order phase transition. 

Let us come back to the questions about the relevance of the unphysical quark degrees 
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of freedom. As already mentioned, deviations from the pure pion gas result become visible 
at T ~ 100 MeV, which corresponds to about 0.6Tc. At this temperature one would not 
expect quark effects to be present in nature. Furthermore, according to universality 
arguments, it is generally believed, that the finite-temperature chiral phase transition in 
QCD with two massless quarks is of second order [Q. This is based on the assumption 
that at Tc there are four massless bosonic degrees of freedom (three pions and one a) 
which determine the infrared behavior of the system. In this case QCD - but also the 
NJL model - should lie in the same universality class as the 0(4)-model, which is known 
to have a second-order phase transition. Although some time ago it was claimed, that 
this argument might not hold if the boson fields are not elementary but composite |^ 



it is probably more likely that the first-order phase transition we observe is an artifact 
of the approximation scheme. In this context the application of renormalization group 
techniques to the NJL model would be extremely interesting. 



7 Conclusions 

We have investigated quark and meson properties within the Nambu-Jona-Lasinio model, 
including meson-loop corrections. These have been generated in two different ways. The 
first method is a systematic expansion of the self-energies in powers ofl/Nc up to next-to- 
leading order [0, In the second scheme, a local correction term to the standard 



Hartree self-energy is self-consistently included in the gap equation |]T4| . We therefore call 
it the "local selfconsistent scheme" (LSS). This scheme can also be derived as the one- 
meson-loop approximation to the effective action Both schemes, the 1/iVc-expansion 



scheme and the LSS, are consistent with chiral symmetry, leading to massless pions in 
the chiral limit. For non- vanishing current quark masses the pion mass is consistent with 
the Gell-Mann-Oakes-Renner relation in the LSS. This is also true in the 1/iVc-expansion 
scheme if one carefully expands both sides of the relation up to next-to-leading order in 

The relative importance of the mesonic fluctuations is controlled by a parameter Am, 
which cuts off the three-momenta of the meson loops. In both schemes we encounter 
instabilities in the pion propagator if the meson effects become too strong. In order to 



find out whether these instabilities are related to an unstable ground state |]19|, leading 
to a "chiral restoration phase transition" at some critical value of A^, we calculated the 
effective action of the LSS for increasing values of Aj\/. (Note, that such investigations are 
not possible within the l/A'^c-expansion scheme, where mesonic fluctuations are included 
only perturbatively.) It turned out, that up to a certain value of Am the system indeed 
seems to move towards a "phase transition" . However, when Am is further increased the 
non-trivial ground state becomes again more stable and no phase transition takes place. 

Of course, at the end, the value of Am, together with the other parameters, has to 
be determined by fitting physical observables. The p-meson and related quantities are 
very well suited for this purpose, since the meson loops are absolutely crucial in order 
to include the dominant p — ^ TTTT-decay channel, while the Hartree-|-RPA approximation 
contains only unphysical gg-decay channels. Here another problem, which constraints 
the possible choice of parameter values, becomes obvious: A priori it is not clear to what 
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extent these unphysical decay modes, which are an unavoidable consequence of the missing 
confinement mechanism in the NJL model, are still present in the region of the p-meson 
peak. 

For the 1/A^c-expansion scheme, the parameters have already been fixed in Ref. ||2CI|| . 
We obtained a reasonable fit of /tt, ('?/''?/') and the pion electromagnetic form factor with 
a constituent quark mass of m = 446 MeV. This means, the unphysical gg-decay channel 
opens at 892 MeV, about 120 MeV above the maximum of the p-meson peak. Further- 
more, the parameters of that fit are far away from the region, where the instabilities in the 
pion propagator emerge. In fact, we found only moderate changes in the pion and quark 
sector: and {il)ip) are lowered by about 20% by the meson loop corrections, while the 
pion mass is increased by about 10%. This indicates that the l/N^. expansion converges 
rapidly and higher-order terms in the l/Ai"c-expansion are small. 

Unfortunately we did not succeed to obtain a similar fit within the LSS. Since in 
this scheme the meson-loop effects lower the constituent quark mass as compared to 
the Hartree mass, it is much more difficult to evade the problem of unphysical gg-decay 
channels in the vicinity of the p-meson peak. We found that a relatively large meson cutoff. 
Km ~ 700 MeV is needed in order to get the quark mass large enough and at the same 
time a fit for However, to our surprise for this cutoff the p-meson self-energy already 
suffers from stability problems, similar to those already discussed for the pion. As a result 
we are not able to get a reasonable description of the p-meson propagator and hence of 
the pion electromagnetic form factor within the LSS. It remains to be checked, whether 
these problems can be cured by taking into account additional intermediate states, like 
vector mesons and axial vector mesons or by different way of regularization. 

In the last part of this article we have investigated the temperature dependence of 
the quark condensate. In both schemes the low-temperature behavior is consistent with 
lowest-order chiral perturbation theory, i.e. the temperature dependence arising from 
a free pion gas. This is a considerable improvement over the mean- field result, where 
the temperature dependence is entirely due to thermally excited quarks, i.e. unphysical 
degrees of freedom. At higher temperatures, however, thermal quark effects also become 
visible in the two extended schemes. We argued that this could be tolerable near the chiral 
phase boundary which is, according to lattice results, identical to the deconfinement phase 
boundary at non-zero temperatures. 

Whereas the perturbative treatment of the mesonic fiuctuations within the l/N^- 
expansion scheme does not allow an examination of the chiral phase transition, this is 
possible in the LSS. For our model parameter set we found a critical temperature of 
164.5 MeV. On the other hand, quark effects are visible already at a temperature of 
~ 100 MeV. Obviously this is still too early to be realistic. Maybe here the model can be 
improved by including additional intermediate meson states. 

In agreement with Ref. |23[ we found a first-order phase transition in that scheme. 
This contradicts the general belief that the non-zero temperature chiral phase transition 
in a model with two light fiavors should be of second order and is probably an artifact 
of the approximation. Here further investigations, e.g. applying renormalization group 
techniques, would be very interesting. 



43 



Acknowledgments 



We are indebted to G.J. van Oldenborgh for his assistance in questions related to his 
program package FF (see http ://www. xs4 al.1.- nl/^ gjvo /i< 'i< '. htrrii ) , which was used in parts 



of our numerical calculations. We also thank G. Ripka, B.-J. Schaefer and M. Urban for 
illuminating discussions. This work was supported in part by the BMBF and NSF grant 
NSF-PHY98-00978. 



A Definition of elementary integrals 



It is possible to reduce the expressions for the quark loops to some elementary inte- 
grals see App. |B| and |C[ In this section we give the definitions of these integrals. 



h = 

Hp)- 

Kip) 

L{PUP2,P3) 



(27r)4 A;2 - m2 + ie 



1 



(27r)4 (F - m2 + te){{k + pf - m? + le) ' 

d'^k 1 

(27r)4 (P - m2 + ieY{{k + pf - + ie) ' 
d^k 1 

(27r)'^ (A;2 — m? + ie){k'l — m? + ie){k2 — m? + ie) 
d^k 1 

(27r)'' [k"^ — m? + ie){kl — m? + ie){k2 — m? + is){k1 — m? + ie) 
d^k kf" 



p^Mi {pi , P2 ) + P2M1 {P2 , Pi 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
,(A.5) 
K,A.6) 



(27r)^ (fc2 — m2 + ie){kl — im? + 'i£){k2 — + i^) 
with ki = k + pi. The function Mi{pi,p2) can be expressed in terms of the other integrals: 
PrP2 Hpi) -pI HP2) + (pI -PrP2) Hpi - P2) + pI {.pl-PrP2) M{pi,p2) 



Mi{puP2) 



2 ((Pi-P2)^ - pIpI) 



(A.7) 



All integrals in Eqs. ( [A.lD to ( [A. 61 ), are understood to be regularized. As described in 
Sec. |5.1| we use Pauli-Villars regularization with two regulators, i.e. we replace 



d^k 
(2^ 



f{k;m) 



(A.8) 



with 



/ij = + J A2 ; Co = 1, ci = -2, C2 = 1 . (A.9) 

One then gets the following relatively simple analytic expressions for the integrals Ji, I{p) 
and K{p): 



(A.IO) 
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(A.12) 



16^ E 2p2(^^.^^^^.^) ( - l-(^.i) - l-(-^.i) + + (A.13) 



with 



4^ 

p2 



(A.14) 



An analytic expression for the three-point function (Eq. |A.4| ) can be found in Refs. [0 
and In certain kinematical regions the four-point function (eq. |A.5| ) is also known 
analytically E^, ESl . 



B RPA propagators 

Using the definitions given in the previous section the gap equation (Eq. ( p.2|) ) takes the 
form 

m = mo + 2igs AN^Nj m h . (B.l) 

Similarly one can evaluate the quark- ant iquark polarization diagrams (Eq. ( |2.6|) ) and 
calculate the RPA meson propagators. The results for a-meson and pion read: 

^"^^^ ^ 1 - 2ig, 2N,Nf (2/i - [p^ - 4m^) I{p)) ' ^^'^^ 

^^^^^ ^ l-2^gs 2N,Nf{2h-p^ /(p)) " ^^'^^ 

If we evaluate these propagators with the constituent quark mass in Hartree approximation 
we can simplify the above expressions with the help of the gap equation (Eq. |B . 1| ) to 
obtain: 

D (v) = (BA) 

"^^^ ^ + 2tgs 2N,Nf (p2 _ 4m2) I{p) ' ^ ' ' 

^^^^^ = ^T2i^:2KN;^jip) ■ ^^-^^ 



As discussed in Sec. this form is also used for the internal meson propagators in the 
LSS. 

A straight-forward evaluation of the vector and axial vector polarization diagrams 
gives 

U,{p) = -t^N.Nf {-2h + {p" + 2m^) I{p)) , (B.6) 
n.,(p) = -t^N.Nf {-2h + (/ - W) I{p)) . (B.7) 
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Because of vector current conservation Up should vanish for p"^ = 0. This is only true if 

m^I{0) = h , (B.8) 

which is not the case if we regularize I{p) and Ii as described in App. 0. This corresponds 
to the standard form of Pauli-Villars regularization in the NJL model 0. Alternatively 
one could perform the replacement Eq. ( ^.11) for the entire polarization loop. In fact, 
this is more in the original sense of Pauli-Villars regularization Then the factor 



in Eq. ( [B.6| ) should be replaced by a factor /i| inside the sum over regulators and one 
can easily show that Eq. ( p.8|) holds (see Eqs. ( |A.1(]| ) and (|A.12|) ). However, this scheme 
would lead to even more severe problems: From the gap equation (Eq. |B.1D we conclude 



that ill should be positive. On the other hand the pion decay constant in the chiral limit 
and in leading order in 1/Nc is given by |^ 

= -2iN,Nf /(O) . (B.9) 

which implies that zm^J(O) should be negative. So irrespective of the regularization 
scheme Eq. ([B.8|) cannot be fulfilled if we want to get reasonable results for m and fi^^ 



at the same time. Therefore we choose the standard form of Pauli-Villars regularization 
in the NJL-model |0 and replace the term Ji in Eq. ( [B.6| ) by hand by -^(0). For 
consistency the Oi is treated in the analogous way. This leads to the following p- and 
Ci-meson propagator 



1 + 2tg,, ^NcNf (-2m2 /(O) + (p^ + 2m2) I{p)) 



^"'^^^ 1 + 2ig^ iNcNf (-2m2 J(0) + (p2 - Am^) I(p)) ' ^^'^^^ 

C Explicit expressions for the meson-meson vertices 

In this section we list the explicit formulae for the meson-meson vertices. We restrict 
ourselves to those combinations which are needed for the calculations presented in this 
article. 

We begin with the three- meson vertices TMi,M2,M3il^P) (^^^ ^^S- §)• 
-2r,,,,,(g,p) = t2mN(^I{p') + I{q) + I{p) + {4m^-^{p'^ + p^ + q^))M{p,-q) 

-iTl'^^Jq,p) = i2mN5ai,{l{p')+p-qM{p,-q)Y 

p2 qfiqX _|_ ^2 pfipX _ p.q (^pf^q^ qf^p^ 



h{q,p) = imN[l{q) + I{p) -21{p') + {Ara^ -2p-q- p'^)M{p,-q) 

= ea,c{q^f{q,p)-p^f{p,q)) , 
f{q,p) = N(^-I{q)+p^M{p,-q) + 2p-qMi{q,-p)^ , (C.l) 
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with p' = -p- q and = AN^Nf. 

For the four-meson vertices we only need to consider the special cases needed for the 
diagrams (b) and (c) in Fig. ^: 

= -n[ ~ + + + 4m^ (M(p, q) + M(p, -g)) 

p^q^ ~i 
+2(m^ (4m^ — p^ — q^) ^) — g, p — g) | 

-^r,,,,,,,(g,p,-p) = -7v|/(p + g) + /(0)+4m2(K(p) + ir(g) + 2M(p,-g)) 

+2p-qM{p, -q) - q^K{q) - p^K{p) 

+m^mm^ - 4p2 - 4g2 + ^^)L{p, -q, 0)| 
-zrj^^,_^(g,p,-g) = (5„6Ar|/(p + g) + /(p-g)+/(4m2-g2)2,(p^_g^p_g)| 

-tTl'^^^Jq,p, -p) = SabN{ - lip + q)- /(O) - (4m2 - q'){Kiq) - p^L{p, -q, 0)) 

+p^K{p) - 2p-q M{p, -g)| 
-iTff,^Jq,p,-q) = -NKabcd{l{p + q) + Iip-q)-p\''Lip,-q,p-q)^ 
-^r:^^i,^(g,p,-p) = -NKabcd{l{p + q) + nO)-p'Kip) 

-q^K{q) + 2p-qM{p, -q) + p^q^L{p, -q, 0) | 

-^r^'a.P.al^.P. = -2<5a;,iv| J(p + q) + I{p - g) + 2/(g) - p-g(M(p, -g) - M(p, g)) 

+ (4m2 - 2p2)(M(p, g) + M(p, -g)) 
+m2(8m2 - + 4g2 + ^ ~ ^^'^"^ )L{p, -q,p- g)| 

-zr^^,^^,^(g,p,-g) = -2<5,,iv|-/(p + g)-/(0) + (/-4m2)ir(p) 

+ (g^ + 2m^)K{q) + (4m^ - 2p-q)M{p, -q) 

+m2(8m2 - 2/ + 4g2 - -g,0)| 

-zr^;';J,^^(g,p, -g) = 2NKa,cd{ - /(p + g) - /(p - q) - 21{q) 

+2p\M{p, q) + M(p, -g)) + p-q{M{p, -g) - M(p, g)) 
+ (2m^p^ -p"^ + {p-qY)L{p, -q,p - g)| 

-^r^;';lp(g,P,-g) = 2NK,,^{lip + q) + I{0)-p'K{p)-{q' + 2m')Kiq) 

+2p-qM{p, -q) + p^{2m^ + q^)L{p, -q,0)^ , (C.2) 
with rp,A/,M,p(g,P, -g) = 9t,u^p'^M,MA^'P'~^'>'^P^M,p,Mi(l^P^-p) = 9,,u^p'^M,p,Mil^P^-p) 

and Kabcd = '^afe'^cd + ^ad^bc — ^ac^bd- 
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D Expressions at non-zero temperature 



To determine the temperature dependence of various quantities we need for the calculation 
of the quark condensate at non-zero temperature in Sec. we adopt the imaginary time 



or Matsubara formalism (see e.g. Ref. ||4^). In principle this amounts to replace the 



integration over energy in the zero temperature expressions by a sum over fermionic or 
bosonic Matsubara frequencies uJn'- 

With this replacement prescription we can define the temperature analogue to the ele- 
mentary integrals, e.g. 

, f d^k 1 

Hp) 



(27r)4 (P - m2 + ie){{k + pf - + ie) 

I{iui,p) = iTY ! ^ ^ (D.2) 

^J{2ny{{iuJn)^-k^-m^){{iuj,, + iuJiy-{k + p)^-m^) 

This example also illustrates our notation: At non-zero temperature the integral de- 
pends on energy and three-momentum separately, which is indicated via a second argu- 
ment. In that way it can be clearly distinguished from its vacuum counterpart with only 
one argument. A similar notation is used for other momentum dependent integrals. The 
non- zero-temperature analogue to the integral Ji will be denoted by Iit- 

We will now summarize the explicit expressions for various temperature dependent 
quantities which are related to the determination of the quark condensate at non-zero 
temperature. The temperature analogue to the gap equation Eq. ( [B.l| ) is given by 

mr = mo - 2gAN,NfmrT E / ' _ j,. 

= mo + 2gAN,NfmThT , (D.3) 



with E = y k'^ + rn^ and u;„ = {2n + l)7rT being fermionic Matsubara frequencies. 
The polarization functions for the RPA mesons read 

li^iiujup) = HN.NfliT - 2iN,Nf{{iuJif -f- Aml)I{iojup) 

li^uup) = AiN,NfhT-2iN,Nf{{iuJif -f )l{iujup) , (D.4) 

with uji = 21ttT being bosonic Matsubara frequencies. Below the phase transition the inte- 
gral IiT can again be replaced with the help of the gap equation Eq. ( p.3|) (cf. Eqs. ( p.2|) 
to O). 

Finally, the constant is given by 

d^p 



D^{iu;i,p){2I{zu;i,p) + I{0,0) - {{luoif - p^ - Aml)K{iuJup)) 
+D^{iuJi,p) ( 3/(0,0) - 'i{{iujif-f) K{iuup) )] , (D.5) 
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where ooi are again bosonic Matsubara frequencies. 
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